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Chapter  1.  Introduction 


There  are  many  practical  problems  that  require  the  understanding  of  turbulence 
behavior  and  therefore  that  call  for  an  analytic  method.  Among  these  practical  problems 
can  be  cited  those  related  to  prediction  of  turbulent  ship-wave  signatures  and  wall 
turbulence  behavior  since  neap  the- wall-effects  are  of  great  importance  for  Reynolds 
average  turbulence  models  for  boundary  layer  models.  There  are  also  problems  related  to 
turbulence  convected  through  the  propulsor  of  a  submarine  or  a  rotor-stator  system  where 
the  turbulent  flow  is  distorted  by  many  different  factors  such  as  the  rapid  curvature,  the 
accelerations  and  the  non-penetration  condition  imposed  by  the  blades. 

In  the  particular  case  of  rotor-stator  system  the  incoming  turbulence  is  heavily 
modified  by  the  distortion  introduced  by  the  first  blade  row,  this  distorted  flow  field  is 
combined  with  the  blade  wakes  generating  a  very  complex  and  new  turbulent  field,  which 
will  interact  with  the  secondblade  row  generating  noise  and  vibrations  (see  figure  1.1). 

In  order  to  predict  and  minimize  noise  and  vibrations  it  is  necessary  to  estimate 
the  two-point  space-time  correlation  function  of  the  turbulence  as  it  goes  through  the 
blades.  To  get  the  two -point  space-time  correlation  function  both  conventional  CFD 
(Computational  Fluids  Dynamics)  and  LES  (Large  Eddy  Simulation)  are  of  little  use. 
CFD  simply  does  not  provide  this  kind  of  correlation  information,  and  LES  is  impractical 
for  realistic  configurations,  which  are  very  complex  and  of  too  large  Reynolds  number  to 
be  modeled  in  this  way.  Perhaps  the  best  way  available  to  obtain  the  twopoint  space- 
time  correlation  function  is  Rapid  Distortion  Theory  (RDT). 

Rapid  Distortion  Theory  tries  to  predict  the  effects  of  a  distortion  imposed  on  a 
turbulent  flow.  When  a  fluid  particle  suddenly  encounters  new  boundaries  (for  instance, 
when  entering  a  contraction  section  of  a  wind  tunnel  or  rotor  passages),  the  sudden 
introduction  of  new  boundaries  modifies  the  mean  flow  field  generating  strains  that  will 
alter  the  turbulence.  When  the  distortion  process  takes  place  in  a  very  short  time,  then  the 
interaction  of  turbulence  with  itself  and  viscous  effects  can  be  neglected.  This  implies,  a 
major  simplification  to  the  calculations,  since  the  formula  can  be  linearized.  It  is  in  this 
sense  that  the  distortion  should  be  rapid  enough. 
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1.1  Previous  Work 

Previous  work  in  this  area  can  be  divided  into  experimental  work  and  aialytical 

work. 

1.1.1  Experimental  Work 

There  is  a  lot  of  experimental  work  on  the  influence  of  ffeestream  turbulence  on 
turbomachinery.  Especially  on  turbulence  effects  on  heat  transfer  and  boundary  layer 
transition  in  turbines  and  turbine  cascades  because  of  its  relevance  to  flows  downstream 
of  combustors,  see,  for  instance,  Hoheisel  et  al.  (1987),  Knshnamoorthy  and  Sukhame 
(1989),  Zhang  and  Han  (1994)  and  Hoffs  et  al.  (1996). 

However  there  are  only  a  few  references  to  experimental  work  involving  the 
influence  of  a  blade-row  configuration  on  turbulent  flow  convected  through  it.  This 
means  there  is  a  little  work  related  to  the  evolution  of  turbulence  as  it  passes  through  a 
propulsor  configuration,  For  instance,  Hobson  and  Shreeve  (1993)  and  Hobson  (1996) 
present  measurements  of  Reynolds  stresses  and  length  scales  (from  time  spectra)  in 
turbulence  approaching  a  cascade.  Ames  and  Plesniak  (1995),  present  measurements  of 
turbulence  stresses  and  timescales  downstream  of  a  turbine  cascade  with  inflow 
turbulence.  There  are  also  Gregory,  Smith  and  Cleak  (1992),  Wunderwald  and  Fottner 
(1996),  Bangert  et  al.  (1997)  and  Ames  and  Plesniak  (1997). 

Bangert  et  al.  and  Ames  and  Plesniak  documerted  the  evolution  of  ffee-stream  as 
it  convects  through  a  blade  row.  Their  work  was  oriented  to  non-homogeneous 
turbulence  evolving  through  a  cascade  (Ames  and  Plesniak)  and  to  rapidly  decaying  free 
stream  turbulence  away  from  the  wall  (Bangert).  Since  thar  work  dealt  with  a  different 
kind  of  flow  than  the  one  for  the  present  study,  it  is  of  limited  relevance  here.  However, 
they  do  show  some  RDT-type  effects  such  as  the  reduction  in  stream-wise  turbulence 
intensity  and  increase  in  length  scales  as  the  flowgoes  through  the  cascade  and  the 
reduction  in  intensity  of  the  wall-normal  turbulent  velocity  component  in  the  vicinity  of 

the  blades. 

If  analytic  studies  of  turbulent  flow  distorted  by  a  cascade  as  the  flow  convects 
through  it  are  to  be  performed,  some  detailed  experimental  data  is  needed.  These 
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experimental  data  should  cover  all  the  elements  affecting  turbulence  evolution  through 
the  cascade  such  as  blade  blocking  effects,  tip  leakage  vortices,  blade  wakes  and 
boundary. 

1.1.2  Analytical  Work 

Many  methods  have  been  developed  for  theoretical  prediction  of  propulsor  effects 
on  turbulence.  Fang  and  Atassi  (1993),  Scott  and  Atassi  (1995),  Silowsky  and  Hall 
(1996),  Peake  and  Kershen  (1997),  present  prediction  methods  for  unsteady  loading  and 
sound  radiation  produced  by  a  blade  row,  including  the  distortion  of  the  turbulence  by  the 
mean  flow  as  it  passes  through  the  cascade.  Majumdar  and  Peake  (1998),  worked  on  the 
distortion  of  turbulence  by  streamline  contraction  that  occurs  upstream  of  a  fan/propulsor 
system.  Graham  (1998),  worked  on  the  distortion  that  turbulent  flow  suffers  as  it  passes 
through  a  cascade  of  flat  plates  at  zero  angle  of  attack  as  a  consequence  of  the  non¬ 
penetration  condition  (or  blade  blocking  effect)  imposed  at  the  blade  surfaces. 

All  this  work  is  ultimately  based  on  Batchelor  and  Proudman's  (1954)  RDT  work. 
Batchelor  and  Proudman  applied  the  linearized  version  of  Cauchy's  equations  to  predict 
the  velocity  wave-number  spectrum  of  turbulence  suffering  homogeneous  distortion 
given  the  RDT  assumptions.  Hunt  (1973),  extended  this  model  to  include  non 
homogeneous  incoming  turbulence  and  nonpenetration  condition. 

In  RDT  the  non-penetration  condition  at  a  solid  surface  may  be  introduced  in 
terms  of  an  irrotational  field  that  modifies  turbulence  intensities  and  it  is  related  to  the 
unsteady  pressure  field.  This  is  why  RDT  has  been  developed  for  aero/hydreacoustic 
predictions  for  rotor  or  stator  noise  generation.  Golubev  and  Atassi  (1993),  Atassi  et  al. 
(1993),  Scott  and  Atassi  (1995),  Silowski  and  Hall  (1998),  Lorence  and  Hall  (1996),  Peak 
and  Kerschen  (1997),  among  others  have  worked  on  this  problem. 

However,  things  are  quite  different  when  trying  to  solve  for  already  disOrted 
turbulence  impinging  on  a  second  set  of  blades,  such  as  a  rotor-stator  system.  In  cases 
like  this  it  is  necessary  to  predict  the  evolution  of  turbulence  as  it  is  distorted  by  the  first 
blade  row.  RDT  has  barely  been  developed  for  treating  such  problems.  Among  the 
researchers  that  have  worked  on  it  can  be  cited  Graham  (1998),  following  Hunt  and 
Graham  (1978)  and  Kullar  and  Graham  (1986)  who  have  worked  on  a  non-staggered 
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cascade  of  flat  plates  at  zero  angle  of  attack.  In  Graham's  (1978)  experiment,  the  blades 
were  separated  by  a  small  distance,  so  that  all  the  distortion  suffered  by  the  homogeneous 
incoming  turbulence  was  due  to  blade  blocking.  He  found  that  the  normalto -the- wall 
component  of  the  turbulent  velocity  field  was  well  predicted  by  RDT,  wheeas  the 
tangential  component  was  not  predicted  that  well.  He  attributed  these  discrepancies  to  the 
pressure  fluctuations  around  the  leading  edge  of  the  blades,  which  were  not  modeled,  and 
to  some  non-linear  evolution  of  the  turbulent  flow  in  the  experiment 

It  could  be  mentioned,  also,  Uzkan  and  Reynolds  (1967),  Thomas  and  Hanckok 
(1977),  Perot  and  Moin  (1995),  and  Walker  et  al.  (1996),  who  have  been  researching 
mainly  on  the  blocking  introduced  by  a  boundary  on  the  free  stream.  Perot  and  Moin 
showed  that  RDT  could  only  model  the  initial  influence  of  the  boundary  upon  the  flow. 
In  the  case  of  a  cascade  configuration  this  means  that  RDT  could  only  be  applied  to 
predict  the  blade  blocking  effect  near  the  leading  edge  region.  However,  many  of  these 
studies  show  that  RDT  provides  good  predictions  for  longer  distances.  As  discussed  by 
Kevlahan  and  Hunt  (1997),  this  could  be  the  result  of  the  cancellation  of  the  nonlinear 
term  specially  when  weak  straining  are  present. 

1.2  Aim  of  the  Present  Work 

Bearing  in  mind  that  this  work  represents  the  first  step  of  a  major  project,  whose 
goal  is  to  provided  a  method  for  noise  prediction  in  a  rotor-stator  system,  recalling,  as 
well,  that  RDT  seems  to  be  the  only  viable  way  to  achieve  this  goal  and  that  there  is 
almost  no  work  (neither  experimental  nor  analytical)  on  turbulence  convected  through  a 
propulsor  type  configuration,  the  main  objective  of  this  work  is: 

•  To  experimentally  characterize  the  evolution  of  a  turbulent  flow  as  it 
convects  through  an  idealized  highly  staggered  cascade  formed  by  non 
symmetric  airfoils. 

•  To  evaluate  the  suitability  of  RDT  in  predicting  the  evolution  of  the  same 
flow  and,  in  so  doing,  to  set  the  basis  for  improvements  in  the  theory. 

To  achieve  this  goal,  measurements  describing  turbulence  convected  through  a 
cascade  were  taken,  RDT  predictions  of  the  same  flow  were  performed  and,  finally,  a 
comparison  between  measurements  and  predictions  was  made. 
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The  Virginia  Tech  Cascade  wind  tunnel  was  used.  Detailed  measurements  of  the 
velocity  components  of  the  flow  field  were  taken  by  using  hotwire  anemometry. 

To  take  the  measurements  it  was  used: 

•  The  Virginia  Tech  Cascade  Wind  Tunnel 

•  Four  Sensors  Hot-Wire  Anemometers,  a  Thermocouple  &  a  Pitot  Static 
Probe 

•  A  Traverse  System  &  Probe  Holders 

•  A  Grid 

2.1  Virginia  Tech  Cascade  Wind  Tunnel 

Figure  2.1  shows  a  general  view  of  the  cascade  wind  tunnel.  The  tunnel  was 
powered  by  a  15h.p.  motor  with  a  fan;  this  is  the  blower  section,  which  was  1.12m  in 
diameter.  Next  to  the  blower  section  there  is  the  expansion  and  a  settling  chamber 
section.  The  total  length  of  both  sections  is  4.26m.  At  the  end  of  the  settling  chamber 
there  are  conditioning  screens,  which  were  intended  to  reduce  turbulence.  Then  proceed 
the  contraction  section  and  the  test  section,  the  last  one  divided  into  two  parts,  the  inlet 
section  and  the  downstream  section,  having  a  total  length  of  5.01m. 

2.1.1  Test  Section 

In  figure  2.2  a  detailed  schematic  of  the  test  section  is  shown.  As  it  was 
mentioned  the  test  section  is  divided  into  two  parts,  the  inlet  and  the  downstream  section. 

The  first  part  has  a  rectangular  cross  section  of  0.76m  x  0.3048m.  This  gives  a 
contraction  ratio  at  the  entrance  of  the  test  section  of  4.8:1.  It  has  a  steel  frame  made  from 
c-section  bars.  This  frame  was  bolted  to  the  floor  to  reduce  vibration  of  the  structure.  The 
lower  end-wall  (see  picture  2.3)  was  made  from  0.75”  of  fin-form  plywood  with  a  smooth 
epoxy  surface  finish  on  top  of  it. 

Right  at  the  beginning  of  the  inlet  section  and  on  the  lower  endwall  were  located 
three  static  pressure  ports.  After  these  pressure  ports  proceeded  a  bleeding  section 
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(described  in  section  2. 1.6)  placed  at  some  0.53m  from  the  contraction,  and  a  suction  slot 
(suction  slots  will  be  described  in  section  2.1.5). 

Both  sidewalls  were  made  from  0.375”  thick  plexiglass.  Both  sidewalls  had  two 
of  these  sheets;  the  first  one  was  screwed  to  the  frame  whereas  the  second  one  was 
attached  to  the  first.  This  way  the  rooftop  had  a  step  to  lie  on. 

The  rooftop  (see  figure  2.4)  was  divided  into  two  parts  a  rectangular  one  and 
triangular  one;  both  of  them  reinforced  with  osection  aluminum  bars.  A  small  slot  of 
around  0.1m  long,  necessary  for  the  reference  static  pitot  probe,  was  cut  on  the  triangular 
part  of  the  rooftop.  Another  slot  of  0.762m  x  0.022m  was  opened  at  the  entrance  of  the 
test  section  where  the  grid  slid  through. 

The  downstream  section  has  a  rectangular  cross  section  of  1.62m  x  0.254m.  For 
the  lower  end  wall  a  0.0625”  Teflon  sheet,  over  a  10/1 6”  of  fin-form  plywood,  was  used 
underneath  the  blade  row.  0.75"  fin-form  plywood  was  used  downstream  the  Teflon 
sheet.  As  before  the  sidewalls  were  made  from  0.375"  plexiglass  butunlike  the  inlet 
section  they  were  single  sheeted,  and  they  are  adjustable  in  order  to  help  in  the  tunnel  set 
up.  Once  the  necessary  orientation  for  the  set  up  was  reached  a  couple  of  clamps  were 
used  to  hold  them  in  place.  The  rooftop  had  two  parts.  One  of  these  parts  has  a  series  of 
slots  drilled  on  it  but  it  is  identical  in  dimensions  to  the  second  and  therefore 
interchangeable  as  required  by  the  measurements.  Gsection  aluminum  bars  reinforced 
both  these  parts.  Every  section  of  the  rooftop  at  the  inlet  and  downstream  sections  were 
held  in  place  by  weights.  See  figure  2.5. 

Finally  three  screens  were  located  at  the  back  of  the  tunnel,  these  three  screens 
were  intended  to  set  up  the  backpressure  for  correct  operation  of  the  suction  slots. 
Aluminum  frames  aid  screens  materials  are  used  for  these  exit  plane  screens.  In  addition 
and  as  a  final  tuning,  two  1"  wide  strips  of  1.62m  long  and  three  1”  vertical  stripes  were 
attached  equally  spaced  along  the  screens  (see  figure  2.6). 

Three  more  static  pressure  ports  were  embedded  on  the  lower  end  wall  at  the  exit 
plane.  Scotch  tape  was  used  around  the  ports  edge  to  ensure  smooth  flow  over  them. 
Figure  2.7  shows  schematics  of  one  these  pressure  ports. 

2.1.2  Blade  Row 
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The  blade  row  leading  edge  was  positioned  at  0.44m  downstream  (in  the  flow 
direction)  of  the  suction  slot  leading  edge.  It  consists  of  8-cantilevered  G.E.  rotor  B 
section  blades  defining  7  passages.  The  blade  maximum  thickness  ratio  of  4.3%  is 
located  at  60%  chord.  Both  the  leading  edge  and  the  trailing  edge  are  rounded.  Figure  2.8 
shows  the  blade  section. 

The  blades  have  a  chord-length  of  0.254m  and  an  aspect  ratio  of  1.  They  were 
tripped  to  ensure  the  flow  underwent  transition  to  turbulence  keeping  the  flow  from 
detaching.  The  trip  strips  were  located  at  0.0254m  from  the  leading  edge  of  each  blade 
and  they  run  all  along  the  span,  they  were  made  from  0.25”  strips  of  0.02”  diameter  glass 
beads.  Every  blade  is  screwed  to  a  0.4m  x  0.1m  aluminum  box  sectioned  support 
structure,  which  allows  the  blades  to  be  moved  at  the  same  time.  In  this  case  it  was  not 
necessary  to  reposition  the  blades.  The  number  of  blades,  according  to  Moore  at  al. 
(1996),  was  enough  to  ensure  a  minimized  effect  from  the  sidewalls  on  the  middle 
passage  (were  the  measurements  were  to  be  performed). 

The  stagger  angle  was  fixed  and  set  at  56.9°,  the  blade  spacing  was  0.236m  and 
the  tip  gap  (gap  between  the  blades  and  the  lower  end-wall)  was  fixed  as  well,  and  set  at 
0.004191m,  the  inlet  angle  of  the  flow  was  65.1°. 

Figure  2.9  shows  a  close-up  picture  of  the  cascade.  There,  passage  number  4, 
between  blades  4  and  5  can  be  seen.  Table  2. 1  shows  the  coordinates  of  a  blade. 


Table  2.1-  Blade  Coordinates  normalized  on  blade  chord  (c=0.254m) 


Lower  Surface 

Upper  Surface 

Tj/c 

& 

T\/c 

0.000000 

0.000000 

0.000000 

0.000000 

0.000435 

0.000596 

0.000060 

-0.001491 

0.001413 

0.001047 

0.000923 

-0.003169 

0.002926 

0.001323 

0.002598 

-0.005009 

0.004966 

0.001388 

0.005091 

-0.006975 

0.007524 

0.001209 

0.008414 

-0.009021 

0.010599 

0.000777 

0.012579 

-0.011102 
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0.014200 

0.000137 

0.017595 

0.019048 

-0.000748 

0.023465 

0.029117 

-0.002550 

0.030187 

0.039178 

-0.004300 

0.037745 

0.049233 

-0.006001 

0.054855 

0.096961 

-0.013419 

0.093151 

0.144562 

-0.019783 

0.140592 

0.192069 

-0.025156 

0.188155 

0.239468 

-0.029599 

0.235822 

0.286809 

-0.033171 

0.283572 

0.334100 

-0.035929 

0.331389 

0.381356 

-0.037929 

0.379254 

0.428588 

-0.039220 

0.427156 

0.475794 

-0.039826 

0.475098 

0.522983 

-0.039750 

0.523069 

0.570167 

-0.038991 

0.571058 

0.617353 

-0.037568 

0.619059 

0.664516 

-0.035603 

0.667097 

0.711679 

-0.032997 

0.715151 

0.758887 

-0.029596 

0.763179 

0.806192 

-0.025241 

0.811130 

0.853654 

-0.019769 

0.858947 

0.901342 

-0.013007 

0.906564 

0.949328 

-0.004778 

0.953911 

0.959464 

-0.002843 

0.963827 

0.969617 

-0.000834 

0.973727 

0.979787 

0.001253 

0.983610 

0.989977 

0.003419 

0.993477 

0.993047 

0.004088 

0.996438 

0.997043 

0.003561 

0.999467 

-0.013180 

-0.015238 

-0.017291 

-0.019400 

-0.021590 

-0.033478 

-0.043940 

-0.053027 

-0.060789 

-0.067278 

-0.072544 

-0.076640 

-0.079613 

-0.081487 

-0.082262 

-0.081938 

-0.080492 

-0.077670 

-0.073277 

-0.067158 

-0.059163 

-0.049143 

-0.036954 

-0.022461 

-0.019107 

-0.015645 

-0.012072 

-0.008389 

-0.007260 

-0.004667 
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1.000000 

0.000000 

1.000000 

0.000000 

2.1.3  Exit  Plane  Screens 

Figure  2.6  shows  the  three  screens  before  mentioned  placed  at  the  ait  plane  of 
the  cascade  located  at  some  1.78m  from  the  blade  row  trailing  edge  (see  picture  2.1). 
These  screens  were  constructed  using  an  aluminum  frame  and  a  screen  material  of  an 
open  area  ratio  of  69.5%.  The  screens  were  clamped  to  the  frame  of  the  tunnel.  They 
were  used  during  the  calibration  process  to  adjust  the  backpressure  to  ensure  the  proper 
working  of  the  suction  slots.  Two  horizontal  and  three  vertical  1"  masking  tape  strips 
were  added  as  a  final  tuning. 

2.1.4  Pressure  Ports 

As  mentioned  before  aluminum  pressure  ports  were  embedded  at  both  the  inlet 
and  downstream  sections  along  the  cross  section,  those  were  used  to  calibrate  the  tunnel 
for  the  first  time  but  they  were  not  used  for  the  present  work.  These  pressure  ports  were 
screwed  to  the  lower  end-wall  and  the  edges  scotch  taped  to  ensure  smooth  flow  over 
them.  The  connection  between  the  ports  and  the  pressure  transducers  was  via  1/16” 
diameter  plastic  tubing.  A  total  of  nine  tabs  were  used;  three  at  the  entrance  of  the  inlet 
section,  three  before  the  suction  slots,  and  three  at  the  exit  plane.  Picture  2.7  shows 
schematics  of  a  pressure  port  embedded  in  the  finform  plywood  floor  with  the  Tygon 
tubing  attached  to  it. 

2.1.5  Suction  Slots 

The  flow  direction  at  the  blade  row  was  65.1°,  which  means  the  boundary  layer 
will  be  thicter  on  one  side  (non-homogeneous  boundary  layer).  Thus  to  get  a 
homogeneous  boundary  layer  at  the  blade  row  leading  edge  line,  two  suction  slots  were 
used  to  remove  the  non-homogeneous  boundary  layer.  Preliminary  studies  indicated  that 
the  boundary  layer  thickness  ranges  from  0.028m  to  0.037m.  To  estimate  the  boundary 
layer  thickness  the  1/7  velocity  distribution  law  was  used.  This  law  leads  to  the  following 
estimation  formula: 
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(  ..  \1/5 


5  =  0.37jc 


4/5 


u 

\  °°s  J 


where  8  is  the  boundary  layer  thickness,  v  is  the  kinematic  viscosity,  U„s  is  the  free 
stream  velocity  in  the  flow-aligned  direction  and  xp  is  the  distance  between  the  inlet 

section  entrance  and  the  point  where  the  boundary  layer  is  being  estimated  (in  the  flow 
direction  as  well).  Since  the  suction  slots  height  was  0.0254m  and  the  boundary  layer  at 
the  suction  slots  was  higher  than  that,  then  these  devices  were  not  able  to  completely 
remove  the  boundary  layer  by  themselves. 

The  amount  of  air  pulled  was  adjustable  by  varying  the  opening  at  the  exit  of  the 
suction  slots.  Figure  2.10  shows  schematics  of  the  suction  slot.  The  picture  also  shows 
the  rectangular  cross  section  wire  of  0.1"  x  0.1"  used  to  trip  the  lower  end-wall  boundary 
layer  from  the  suction  slot.  The  wire  was  placed  at  0.182m  axially  (direction 
perpendicular  to  the  blade  row  leading  edge  line)  upstream  of  the  leading  edge  line  of  the 
blade  row  and  0.007m  downstream  of  the  leading  edge  of  the  suction  slot. 

Also  used  during  the  tunnel  set  up  were  a  couple  of  sidewall  openings  formed  by 
the  gap  between  the  leading  edge  of  the  first  and  the  last  blades  and  the  inlet  sidewalls. 
These  openings,  of  about  0.254m  height  and  0.03m  deep,  were  adjusted  by  using  an 
aluminum  flat  platetaped  to  the  exterior  of  the  sidewalls  (see  figure  2.2). 


2.1.6  Bleeding 

Since  the  suction  slots  were  not  able  to  completely  remove  the  boundary  layer, 
then  extra  bleeding  was  necessary. 

Figure  2.11  shows  a  porous  strip  spanning  the  lower  wind  tunnel  end- wall,  to 
bleed  8%  of  the  flow  and  in  so  doing  enabling  the  suction  slots  to  work  without  imposing 
any  pressure  gradient  on  the  flow.  This  bleeding  was  placed  at  0.5334m  from  the  exit  of 
the  contraction  section.  The  strip  was  0.0635m  width  and  0.762m  long  and  made  from 
stainless  steel. 
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2.2  Four  Sensor  Hot-Wire  Anemometers,  Thermocouple  &  Pitot  Static 

Figure  2.12  shows  a  four-sensor  probe,  model  AVOP-4-lOO  manufactured  by 
Auspex  Corporation.  This  probe  is  a  miniature  Kovaznay  type  probe  with  four  sensors 
arranged  in  two  orthogonal  X-wire  arrays  on  eight  stainless  steel  prongs  (75  pm  in 
diameter  at  their  tips)  made  from  etched  tungsten  wire  of  5pm  in  diameter  and  0.8  mm 
long  giving  a  length  to  diameter  ratio  of  160.  The  wires  are  inclined  at  a  nominal  45 
angle  to  the  probe  axis.  The  measuremert  volume  for  the  sensor  was  0.5  mm3.  Every 
sensor  of  the  probe  was  connected  to  a  constant  temperature  anemometer  unit  (Dantec 
model  56C17/56C01),  which  was  optimized  to  give  a  frequency  response  greater  then  20 
kHz.  Hot  wire  signals  were  buffered  by  fourx  10  buck-and-gain  amplifiers  containing 
calibrated  RC  filters  to  limit  their  response  to  50  kHz. 

A  Pentium  PC  via  a  Hewlet  Packard  El 406 A  16-channel/ 16-bit  digitizer  recorded 
the  output  voltages  from  the  anemometer  unit.  This  PC  was  also  used  to  control  the 
traverse  system. 

A  thermocouple  was  used  to  measure  the  flow  temperature.  The  copper- 
constantan  thermocouple,  manufactured  by  OMEGA  Instruments,  was  connected  to  an 
OMEGA  Instruments  DP80  series  temperature  transducer. 

The  free  stream  velocity  was  measured  by  using  a  Pitot  static  probe  connected  to 
a  SETRA  model  239-pressure  transducer  with  an  input  range  of  0-5  in  of  water  and  an 
output  of  1  V/in  of  water. 

2.2.1  Hot  wire  Calibration 

The  probe  was  velocity  calibrated  before  and  after  every  set  of  measurements  to 
account  for  possible  changes  in  the  wire  characteristics.  To  do  so  the  probe  was  placed  in 
a  uniform  jet  of  a  TSI  calibrator  and  King's  law  was  used  to  correlate  the  output  voltages 
to  the  cooling  velocities. 

Velocity  components  were  determined  by  a  direct  angle  calibration  where  the 
probe  was  pitched  and  yawed  over  a  range  of  angles  making  with  those  rnlues  look-up 
tables  to  get  the  relationships  between  cooling  velocities  and  flow  angle;  see  Wittmeref 
al.  (1998).  That  study,  where  measurements  in  a  pipe  were  taken,  revealed  that  the  four 
sensor  probe  gives  accurate  measurements  of  velocities  in  the  nor  the  wall  region 
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(within  10%  error  at  0.0025m  off  the  pipe  wall),  which  is  an  advantage  for  the  cascade 
measurements,  where  the  minimum  distance  off  the  wall  was  0.00508m. 

Figure  2.13  shows  a  hot-wire  probe  placed  in  the  calibration  equipment.  There, 
the  probe  can  be  seen  attached  to  the  pitching-yawing  gear  in  front  of  the  jet  of  a  TSI 
calibrator  model  1125. 

2.3  Traverse  System 

The  traverse  system  (see  figure  2.14)  was  a  two  axes  movement  system.  It 
consisted  of  a  couple  of  lead-screw,  accurate  to  0.0025m  per  meter,  mounted  on  a  double 
rail  and  moved  by  a  couple  of  stepper  motors.  The  stepper  motors  were  manufactured  by 
Compumotor  and  the  model  was  S-57-83-MO.  The  stepper  motors  were  controlled  by  a 
PDX13  single-axis  package  mini-step  drive.  The  resolution  of  the  programmable  traverse 
movement  was  2.54e-5m.  The  lead  screws  moved  a  carriage  where  the  probe  holders 
were  attached.  These  carriages  were  hold  to  the  screws  using  anti-backlash  nuts. 

Technolsel  manufactured  all  the  anti-backlash,  the  lead  screws,  the  double  rail 
and  the  carriage  system.  The  whole  traverse  system  was  mounted  on  an  I-section 
aluminum  bar. 

Two  different  probe  holders  were  used  depending  on  the  in-passage  section  it  was 
measured.  Figure  2.15  shows  the  probe  holders. 

2.4  Grid 

A  turbulence-generating  grid  was  especially  designed  for  these  measurements. 
The  turbulence  was  intended  to  be  isotropic  and  homogeneous  and  should  be  fully 
developed  at  the  blade  row.  To  meet  these  requirements  a  grid  was  specially  designed. 
The  grid  (see  figure  2.16)  consists  of  10  vertical  and  4  horizontal  steel  rods  welded 
together.  The  rods  diameter  was  0.0159m;  the  grid  cell  size  was  0.0749m  and  the  open 
area  ratio  62%.  The  grid  was  placed  at  1.6m  upstream  from  the  center  of  the  blade  row 
leading  edge  line  (see  figure  2.17).  This  distance  was  more  than  enough  to  ensure  fully 
developed  turbulence  at  the  middle  passage,  where  measurenents  were  taken.  According 
to  P.E  Roach  (1987),  the  minimum  distance  down  stream  of  the  grid  necessary  to  get 
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fully  developed  turbulence  was  12  times  the  cell  size  (0.9m)  whereas  the  shortest  distance 
in  the  present  work  between  the  blade  row  and  the  grid  was  1 .04m. 

2.9  Cascade  Set-Up 

Figure  2.2  shows  the  coordinate  system  used  during  the  calibration  process.  The 
origin  was  located  on  the  lower  end-wall  at  the  center  of  the  line  that  runs  from  the 
leading  edge  of  blade  number  4  to  the  leading  edge  of  blade  number  5.  Thez-axis 
(tangential  direction)  is  parallel  to  the  leading  edge  line  and  positive  towards  the  short 
sidewall  (see  figure  2.2).  The  x-axis  (axial  direction)  is  perpendicular  to  the  z-axis  and 
positive  downstream;  the  y-axis  completes  the  right  hand  coordinate  system  being 
positive  upward.  Distances  will  be  normalized  either  on  the  axial  chord  £a  =  0.1386m)  or 
on  the  geometric  chord  (c  =  0.254m);  in  either  case  it  will  be  indiated.  This  system  will 
be  called  axial  tangential  and  will  be  used  through  out  the  present  study. 

While  calibrating,  two  pitot-static  probes  connected  to  a  couple  of  SETRA  model 
239  pressure  transducer  with  an  input  range  of  0-5  in  of  water  and  an  output  of  1  V/in  of 
water  were  used.  The  first  probe  measured  the  dynamic  pressure  (reference  pressure, 
Pref.)  and  it  was  placed  at  0.42m  downstream  of  the  lower  end-wall  bleed;  its  coordinates 
in  axial  tangential  coordinate  system  werex=0.306m,  z=0.6m  andy=0.18m.  The  second 
probe  measured  both  the  static  and  the  total  pressure  relative  to  the  atmospheric  pressure. 
This  is  Px  =  -  Pa  and  P2  =  Psd  -  Pa .  Where  P,d  is  the  total  pressure,  Psd  is  the  static 

pressure  and  Pa  is  the  atmospheric  pressure.  This  probe  was  moved  by  the  traverse  along 
a  plane  located  at  x=-0.127m  ahead  of  the  blade  row  leading  edge  line  and  parallel  to  it  at 
y=0.127m. 

A  thermocouple  was  used  to  measure  the  flow  temperature;  The  thermocouple, 
manufactured  by  OMEGA  Instruments,  was  connected  to  a  temperature  transducer  also 
manufactured  by  OMEGA  Instruments  model  DP80.  It  was  placed  at  some  0.3048m  from 
the  beginning  of  the  inlet  section,  0.05m  from  the  longer  sidewall  and  0.05m  from  the 
upper  end- wall.  Its  axial  tangentid  coordinates  were  x=-0.84m,  z=  1.035m  andy=0.23m. 

The  variables  used  during  the  calibration  were,  the  lower  and  upper  endwall 
suction  slots,  the  inlet  section  sidewall  suction  slots  (one  on  each  sidewall),  the 
downstream  section  sidewall  angles  and  tie  backpressure  (which  will  depend  on  the 
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screens  at  the  exit  plane).  Every  other  gap  on  the  inlet  section  was  sealed  including  the 
port  where  the  moving  probe  holder  came  through. 

At  the  moment  of  these  measurements  a  moving  wall  underneath  the  row  blade 
was  being  used  in  another  experiment,  this  moving  wall  needed  two  gaps  of 
1.4mx0.006m,  one  on  every  sidewall.  The  gaps  were  placed  starting  at  0.1m  downstream 
of  the  lower  end-wall  suction  slot  leading  edge  and  continuing  in  the  downstream 
direction.  These  gaps  were  kept  open  at  every  time. 

Another  important  variable  was  the  turning  angle  of  the  flow  after  the  blade  row. 
This  was  set  at  1 1.78°,  see  Wang  (2000),  for  details  on  the  setting  up  procedure. 

The  objective  of  the  calibration  was  to  get  the  tunnel  working  at  the  desired  flow 
conditions,  that  is  no  mean  velocity  variation  along  the  inlet  section,  periodic  static 
pressure  distribution  and  no  pressure  gradients  across  the  test  section  (especially  at 
passage  4,  where  the  measurements  were  to  be  taken). 

After  many  different  tries  the  best  possible  configuration  is  the  following:  LSG 
covered  with  0.1016mx0.03048m  blocker,  SSG  covered  with  0.1 01 6mx0. 02032m 
blocker,  and  top  and  bottom  suction  slots  opened  0.01905m  (3/4”).  Three  screens  with 
two  horizontal  and  three  vertical  1"  strips  of  masking  tape  attached  to  last  screen  were 
used.  The  horizontal  strips  were  1.62m  long  whereas  the  vertical  strips  were  0.254m 
long.  For  this  case  there  was  almost  no  pressure  gradient  (3%  pressure  difference  across 
the  passage)  and  the  mean  velocity  did  not  change  in  the  inlet  section.  Figure  2.18  shows 
pressure  plots  for  the  best  configuration. 
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Measurements  were  taken  to  document  turbulent  flow  as  it  convects  through  a 
linear  compressor  cascade  containing  many  of  the  features  found  in  an  actual  rotor  such 
as  blade  thickness,  flow  curvature,  deceleration  and  tip  leakage.  Measurements  were 
taken  with  and  without  a  turbulence-generating  grid.  The  reference  system  used  while 
taking  the  measurements  was  defined  as  shown  in  figure  3.1;  the  origin  located  at  the 
lower  end-wall  midway  between  the  two  blades  defining  the  central  passage  (passage  4 
between  blades  4  and  5-see  figure  2.2),  the  x  axis  was  perpendicular  to  the  leading  edge 
line  pointing  downstream,  the  z  axis  was  parallel  to  the  leading  edge  line  pointing 
downward  and  the  y  axis  completing  the  system  according  to  the  right  hand  rule  (positive 
coming  out  of  the  lower  end- wall).  This  reference  system  will  be  called  from  now  on  the 
'axial-tangential  reference  system'. 

Data  along  eight  measurement  planes  upstream,  in  and  downstream  of  the  central 
passage  were  gathered,  the  first  two  placed  before  the  leading  edge  line,  five  within  the 
passage  and  the  last  one  after  the  trailing  edge.  As  shown  in  figure  3.1  the  first 
measurement  plane  was  located  at  x=-0.8cfl  ahead  of  the  leading  edge  line  (ca  —0.1387m 
is  the  chord  length  in  the  V  direction).  This  will  be  known  as  'plane  O’.  The  second 
measurement  plane,  'plane  1'  was  located  atx=-0.23  ca  ahead  of  the  leading  edge  line.  The 
third  plane,  'plane  2',  coincides  with  the  leading  edge  line.  Plans  3,  4,  5,  6,  7  were 
located  at  x=0.27ca,  x=0.48cfl,  x=0.77cfl,  *=0.98ca,  x=1.26ca  respectively.  Every  plane 
extended  over  0.127m  in  height  measured  from  the  lower  end  wall  (approximately  half 
the  span  of  the  blades).  This  way  the  flow  before  the  blade  row,  within  the  passage  and 
behind  the  blade  row  was  covered.  A  total  of  about  200  points  per  plane  were  measured. 

The  velocity  components  are  U,  V,  W  and  the  Reynolds  normal  stress  components 
ua,  v'2  and  w*2  are  defined  in  the  x,  y  and  z  directions  respectively.  In  general  distances 
were  normalized  on  the  axial  chord  (ca  =0.1 387m),  mean  velocity  and  turbulent  stress 
components  were  normalized  on  the  approach  free  stream  velocity  which  was  constantly 
monitored  by  a  Pitot  probe  placed  at  0.9271m  downstream  of  the  grid  which  was  placed 
at  1.6m  upstream  from  the  middle  of  the  blade-row  leading  edge  line  in  the  flow 
direction,  0.102m  from  the  upper  end-wall  and  0.35m  from  the  left  side  wall  looking 
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downstream.  In  terms  of  the  axial  tangential  coordinate  system  the  probe  was  placed  at 
x=-0.306m,  z=0.6m  andy=0.18m. 

The  Reynolds  number  based  on  the  blade  chord  of  0.254m  and  on  the  approach 
free  stream  velocity  of  25m/s  on  average  was  370,000;  The  Reynolds  number  based  on 
the  grid  size  (M=0.07493m)  was  1 24,500. 

Uncertainties  of  measured  quantities  computed  at  20:1  odds  are  given  in  table  3.1. 


Table  3.1 


Quantity 

Uncertainty 

U,  V,  w 

±1%  Uoos 

uu 

±1%  u 

vu 

±6%  v} 

w* 

±6%  vv2 

tke 

±3%  tke 

Where  tke  is  the  turbulence  kinetic  energy  computed  as  (w'2+v'2+>v’2)/2  and  Ums  is 
the  approach  free  stream  velocity  in  the  flow  direction.  These  calculations  were  made 
according  to  Kline  and  McClintockfc  (1953)  method.  The  principal  sources  of  error  are 
represented  by: 

•  A/D  converter 

The  random  voltage  uncertainty  in  the  A/D  conversion  was 
3.662E-4  Volts. 

•  Velocity  and  angle  calibration 

Kingfe  law  was  used  to  relate  the  output  voltages  to  the  cooling 
velocities.  The  uncertainty  associated  with  the  curve  fitting  was 
0.013 141 5m/s. 

•  Uncertainties  in  the  reference  velocity 

The  reference  velocity  was  measured  by  means  of  a  Pitot  static 
probe  connected  to  a  pressure  transducer.  The  uncertainties  of  this 
pressure  transducer  were  0.01  in  of  water  or  2.4884Pa. 
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•  Uncertainties  due  to  the  total  number  of  samples  measured. 


3.1  Flow-field  characteristics 

Three  different  regions  can  be  identified:  the  incoming  flow  region  upstrem  of 
the  blade  row,  the  viscous  region  (within  the  passage  and  behind  the  row  blade)  and  the 
two  dimensional  region  where  the  V  component  of  the  mean  velocity  field  in  they 
direction  can  be  neglected. 

The  viscous  region  is  represented  by  the  blade  bounday  layers,  the  lower  and 
upper  end-wall  boundary  layers,  the  tip  leakage  vortex  (where  the  flow  is  three 
dimensional)  and  the  blade  wakes  behind  the  blade  row  (see  figure  3.2). 

In  the  two  dimensional  region  the  mean  velocity  component  in  they  direction  can 
be  neglected  when  compared  to  the  other  two.  This  region  has  been  characterized  through 
the  data  gathered  at  the  first  measurement  plane. 

3.1.1  Incoming  flow 

To  characterize  the  incident  flow,  measurements  atx/ca=-0.8  were  used.  This  way 
the  characteristics  of  the  grid-generated  turbulence,  the  turbulence  of  the  free  stream 
without  the  grid  and  the  boundary  layer  thickness  on  the  lower  end-wall  were  determined. 

Here  for  simplicity,  a  new  reference  system  is  introduced.  The  new  reference 
system,  called  'flow  aligned  reference  system',  is  represented  byr^  in  the  flow  direction, 
zs  perpendicular  to  xs,  parallel  to  the  end-wall  and  pointing  blade  number  four  and, 
finally,  ys  defining  the  right  hand  system  (positive  upward).  Us  is  the  mean  velocity 
component  in  the  xs.  The  other  two  mean  components  are  zero.  The  equations  used  to  get 
the  velocity  components  are: 


u; 

cos6y  COS02 

sin02 

-sinf?^  cos02 

~U~ 

0 

= 

-sin02  cos0^, 

COS02 

-sin02  sin©^ 

V 

0 

sin0. 

0 

COS0, 

w 

where 
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Q 


y 


rWy 

f 

1 

=  -tan  1 

and  Gz  =  tan 

< 

V 

U  cosdy-Wsindy 


A 


J 


the  rotation  sequence  was  dy-62 . 

Figures  3.3a,  b  and  c  show  the  distribution  of  us'/vs'  and  us/ws'  along  the  first 
measurement  plane  between  z=0.3ca  and  z=0.8ca  at  three  different  heights  (0.0762m, 
0.1016m  and  0.127m  from  the  lower  end  wall)  with  the  grid  on.  Here  us\  ws'  and  v/  are 
the  R.M.S.  turbulence  fluctuations  in  xs,  zs  and  ys  directions  respectively.  The  level  of 
isotropy  is  better  than  expected  considering  that  the  grid  was  not  sited  at  the  contraction 
section  of  the  wind  tunnel  and  that  there  was  8%  diffusion  over  the  porous  bleed  slot.  The 
difference  among  turbulent  components  was  no  more  than  5%,  whereas  for  grid 
generated  turbulence  this  percentage  rises  up  to  10%.  The  turbulence  intensity  at  the 
center  of  this  measurement  plane  (z/ca=l  .6)  was  3.1%  of  U<*s.  These  three  plots  also  show 
a  variation  in  the  turbulence  intensity  ratio  across  the  passage.  Since  the  grid  was  placed 
perpendicular  to  the  tunnel  flow  but  the  cascade  and  the  measurement  planes  were  at  a 
large  angle  (65.1°)  with  respect  to  the  vertical  measurements  made  at  differentz  locations 
were  made  also  at  different xs  locations.  Therefore,  this  variation  represents  the  decay  the 
turbulence  suffers  due  to  viscous  dissipatbn.  It  is  interesting  to  note  that  the  decay  is  not 
only  affecting  the  turbulent  velocity  components  but  also  the  ratio  of  them.  This  can  be 
explained  by  the  fact  the  incoming  turbulent  flow  is  not  100%  homogeneous  and 
isotropic.  The  decay  effect  is  also  drown  in  figure  3.3d.  Tke,  averaged  on  values 
measured  at  the  same  three  different  heights  above  mentioned,  has  been  plotted  along  the 
inflow  plane.  The  same  variation  observed  in  the  turbulent  components  is  present  here. 
As  already  explained  this  is  due  to  the  relative  position  between  the  grid  and  the 
measurement  planes.  Thus  this  plot  can  be  used  to  estimate  the  rate  of  dissipation  due  to 
viscosity. 

The  rate  of  dissipation  for  spatially  developing  grid  turbulence  can  be  computed 


as: 


£  =  -1/ 


d(tke) 

dxs 
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where  Us  is  the  mean  velocity  component  in  the  direction  of  convection, 
been  computed  by  a  finite  difference  method  as: 


djtlce) 

dxx 


has 


A  {the)  _  A  (the) 

Axs  Azcos(90-65.1) 


where  A(tke)  is  the  difference  between  two  consecutive  points  on  the  inflow  plane 
separated  by  A z  and  65.1°  is  the  angle  between  the  leading  edge  line  and  the  vertical  (see 
figure  2.2).  Formula  3.1  was  applied  to  every  point  on  the  infbw  plane  at  y=0.5494ca, 
y=0.7326cfl  and  y=0.9157ca  to  get,  then,  e  at  each  one  of  these  points  and  performing, 
afterwards,  an  average  on  all  of  them.  This  calculation  gives  a  dissipation  rate  ofe=- 
0.000198  normalized  on  Uljca.  The  central  postulate  of  the  Kolmogorov  theory  of 

turbulence  states  that  if  the  Reynolds  number  is  sufficiently  high  there  exists  a  scaling 
region,  the  inertial  sub-range,  which  separates  the  low  wave-number  region  from  the  high 
wave-number  region  in  which  viscous  dssipation  occurs.  If  this  sub-range  exists  then 
Euusiki),  which  is  the  one-dimensional  wavenumber  energy  spectrum  in  the  stream-wise 
direction  related  to  u's,  can  be  computed  as 


Euus  =Cle3eki3  (3.2) 

where  C,  is  the  Kolmogorov  constant, 

C,=  4.5 -2.4/7  (3-3) 

and  n,  the  slope  of  the  inertial  sub-range,  is  5/3  for  locally  isotropic  turbulence. 

For  the  present  flow  C,,  estimated  from  measurements  using  equation  3.2,  at 
every  measured  point  at  three  different  heights  (j^O. 0762m,  y^O.  1016m  andy=0.127m) 
and  then  averaged,  is  0.2.  With  this  value  of  C,  and  using  equation  3.3,  n  turned  out  to  be 
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1.79,  which  is  away  from  the  ideal  value  of  5/3.  This  means  that  the  inertial  subrange  is 
smaller  and,  therefore,  the  region  of  the  spectrum  under  the  influence  of  viscous 
dissipation  is  bigger,  which  brings  up  the  idea  that  viscous  effects  could  be  significant  in 
this  flow. 

Consistent  with  the  level  of  homogeneity  and  isotropy  of  this  turbulent  flow, 
figure  3.3e  shows  the  level  of  agreement  between  the  three  components  of  the  spectrum 
computed  with  the  von  Karman  interpolation  formula,  (see  section  4.2,  formula  4.21), 
and  the  three  components  of  the  measured  time  spectrum  (measured  at  y=0.101m, 
2=0. 1345m)  at  frequencies  within  and  below  the  inertial  subrange. 

The  frequency  spectra  are  defined  as 

G.„  =  2ijmI&[x,(/,r)|2] 

where  Ex  is  the  expected  value  of  \Xk(f,T)\  which  is  the  Fast  Fourier  Transform  (FFT) 

of  the  turbulent  velocity  component  u's  and  T  is  the  time  interval.  Similar  expressions 
exist  for  Gws  and  .  In  terms  of  units  the  autospectra  components  are  defined  as 

HI2 

s2Hz 

The  spectra  are  presented  in  flow-aligned  coordinate  system  and  normalized  on 
U2oos.  They  were  measured  at  a  sampling  rate  of  51.2kHz  by  averaging  the  Fourier 
transforms  of  50,  1024-point  records.  The  lowest  resolvable  frequency,  just  a  littleabove 
of  50  Hz,  seems  to  be  in  the  flat  part  of  the  spectrum.  So  this  can  be  used  to  estimate  the 
integral  length  scales.  Taylor's  frozen  hypothesis  states  that  the  time  distribution  of  a 
turbulent  quantity  at  a  point  is  nearly  identical  to  the  instantan©us  distribution  of  the 
same  quantity  along  an  axis  through  that  point  in  the  flow  direction.  According  to  this  the 
stream-wise  integral  length  scales,  which  represents  the  distance  over  which  the  influence 
of  turbulence  extends,  can  be  calculated  as: 
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A«=^£«(| f)=^c“(50) 

a  ^  I?  /27r50^_  n  szrw 

^ vs  *2  ^vws(~~rr  )  “  >  ’2  ' 

2w5  4ws 

a  _  ^  i7  / 27T 50^  __  U s  ^  f^cw 

^ws  ”  r~T ^wwj(  r7  a  '2^whwV  ) 


where  Euus(^^-),  Exvs(^p-)  and  E^&p-)  are  the  corresponding  one  dimensional 

wave-number  spectra  for  50Hz,  Since  the  lowest  frequency  bin  was  used,  this  length 
scales  represent  the  maximum  distance  over  which  the  influence  of  turbulence  extends. 
The  wave  number  spectra  are  related  to  the  time  auto  spectra  as  follows: 


E  =— hhTjl  similar  for  the  other  two.  In  terms  of  units  the  wavenumber  spectra 

u“  2k 

m2 

components  are  defined  as  — r — ,  where  ki=2nf/Us  and  /is  the  frequency  in  Hz. 

s  kx 

For  homogeneous  and  isotropic  turbulence  the  longitudinal  integral  length  scale, 
Am  related  to  u's,  should  double  the  integral  lateral  length  scales,  Aw related  to  v's  and 
A  ^related  to  w's.  This  is  shown  in  figure  3.3f.  Equation  3.4  gives  a  Aus  of  0.022m  which 
is  16%  of  the  axial  chord  length  and  22%  of  the  upstream  projected  blade  spacing, 
computed  as  0.235m/cos(65.1 ). 

Without  the  turbulence-generating  grid  turbulence  levels  are  much  lower  and 
incoming  turbulence  is  not  isotropic  anymore.  After  extracting  contributions  due  to 
electrical  noise,  R.M.S  of  turbulence  fluctuations  us',  w/  and  v/  turned  out  to  be  0.12%, 
0.21%  and  0.23%  of  U„s  whereas  the  corresponding  integral  length  scales  Am ,  Aws  and 

Avs  were  0.00838m,  0.0168m  and  0.0122m.  Predictions  without  the  grid  were  not 

performed  since  the  hi$i  level  of  anisotropy  implies  an  inflow  wavenumber  spectrum 
that  is  difficult  to  define.  A  welldetermined  incoming  flow  in  terms  of  the  turbulence  is 
very  important  since  this  work  represents  the  initial  steps  in  understanding  RDT,  so  the 
flow  has  to  be  kept  simple. 
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Measurements  of  the  boundary  layer  were  performed  but  not  in  detail  since  the  4 
sensor  probe  used  was  not  appropriate  for  this  purpose.  However,  the  objective  here  was 
to  check  the  tangentially  uniformity  of  the  lower  endwall  boundary  layer.  The 
measurements  revealed  that  the  boundary  layer,  with  and  without  the  grid,  was 
tangentially  uniform  and  its  thickness  was  0.01524m  at  this  measurement  plane. 


3.1.2  Viscous  region 

As  mentioned  at  the  beginning  of  section  3.1,  the  viscous  regions  in  a  linear 
compressor  cascade  are  formed  by  the  lower  and  upper  end-wall  boundary  layers,  blade 
boundary  layers,  blade  wakes  and  most  important  by  the  tip  leakage  vortex.  For  he 
present  work  the  description  of  the  viscous  region  and  mainly  of  the  tip  leakage  vortex  is 
important  to  identify  the  zones  where  three-dimensional  effects  dominate  the  flow.  Once 
this  three  dimensional  region  has  been  identified  the  two  dimensional  regin,  where  the 
comparison  between  measurements  and  predictions  will  be  performed,  can  be 
determined.  Two-dimensional  implies  that  only  two  velocity  components,  those  inx  and  z 
directions,  define  the  velocity  field. 

At  this  point,  a  new  reference  system  is  defined.  This  reference  system  is  aligned 
with  the  potential  core  of  the  flow  in  the  two  dimensional  region.  The  two  dimensional 
RANS  solution  (see  section  4. 1.3.2)  of  Shin  et  al.  (1999)  for  the  present  tunnel 
configuration  was  used  to  get  the  streamline  coordinates  through  the  center  of  the 
passage,  necessary  to  determine  the  relationship  between  the  axial  tangential  coordinate 
system  and  the  'center  stream  line  aligned  reference  system'.  The  center  streamline 
aligned  reference  system  is  represented  byxc>  positive  downstream,  zc,  perpendicular  to  xc 
pointing  to  the  pressure  side  of  blade  4  and  yc,  completing  the  system  according  to  the 
right  hand  side  rule. 

The  relationship  between  these  two  reference  systems  is: 


UC=U cos6yc  -  W sin0K 
VC=V 

Wc=Wcos6yc+Usindyc 
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-*(W\  U  and  W 


KUJ 


where  6yc  is  the  rotation  angle  around  yc  and  it  was  computed  as  tan 

are  the  velocity  components  of  the  flow  at  the  center  streamline.  Values  of  dyc  and  the 
coordinates  at  where  these  quantities  were  computed  are  listed  in  table  3.2 


Table  3.2:  Rotation  angles 


Measurement  Plane 

(x/cj 

Center  streamline*  /  ca 

coord. 

Center  streamlinez  /  ca 

coord. 

-0.8 

64.74° 

-0.8 

1.69 

-0.23 

66.13° 

-0.23 

0.45 

0 

62.18° 

0 

0 

0.27 

58.07° 

0.27 

-0.49  ; 

0.48 

55.85° 

0.48 

-0.82 

0.77 

54.18° 

0.77 

-1.24 

0.98 

53.4° 

0.98 

-1.5 

1.26 

52.51° 

1.26 

-1.9 

The  viscous  region  as  well  as  the  influenceof  the  turbulence-generating  grid  can 
be  visualized  through  the  following: 

-)  Mean  stream-wise  velocity  (Us/UooS) 

-)  Maximum  mean  velocity  deficit  with  respect  to  the  local  center  streamline 
aligned  velocity  component  outside  the  tip  leakage  vortex. 

-)  Mass  flow  deficit  with  respect  to  the  potential  core. 

-)  Mean  stream-wise  vorticity 

-)  Turbulence  kinetic  energy  (fke) 

Figure  3.4  a  to  h  shows  contour  plots  of  mean  stream- wise  velocity  normalized  on 
Uoos  at  planes  4,  5,  6  and  7,  with  the  grid  (a,  b,  c  and  d)  and  without  the  grid  (e,  f,  g  and 
h).  As  mentioned  before  around  200  points  has  been  measured  at  every  plane  and  an 
interpolation  has  been  performed  afterwards  in  order  to  generate  these  contours.  Only  the 
last  four  measurement  planes  are  shown  since  nothiig  much  that  could  help  determining 
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the  two-dimensional  region  can  be  seen  in  the  first  three.  In  these  plots  the  observer  is 
positioned  upstream,  looking  downstream.  The  axes  of  the  plot  are  axial  tangential.  This 
method  will  be  applied  on  every  cross  secton  contour  plot  along  the  present  work. 

As  shown  by  this  set  of  plots  the  mean  velocity  field  distribution  is  similar  with 
and  without  the  grid.  However  the  size  of  the  zone  under  the  influence  of  the  tip  vortex  is 
about  18%  in  average  less  with  out  the  grid  than  with  the  grid.  The  area  under  the 
influence  of  the  tip  leakage  vortex  was  bounded  at  its  top  by  a  line  at  which  the  velocity 
was  considered  to  be  99%  of  the  local  velocity  of  the  potential  core.  Under  this  criterion 
the  region  influenced  by  the  ip  leakage  vortex  with  and  without  the  grid  was  estimated 
and  compared  to  each  other. 

A  region  of  large  velocity  deficit  can  be  seen  near  the  lower  endwall  region 
revealing  the  presence  of  the  tip  leakage  vortex.  The  tip  leakage  vortex  is  formed  as  the 
result  of  the  flow  passing  through  the  gap  between  the  blade  and  the  end-wall  due  to  the 
pressure  difference  between  the  suction  and  the  pressure  side  of  the  blade. 

As  the  flow  goes  downstream  the  vortex  increases  in  size  reaching  a  maximum 
height  of  about  0.4ca  at  plane  7,  and  it  moves  to  the  pressure  side  of  the  passage. 

These  pictures  also  show  a  velocity  gradient  across  the  passage  width  with  lower 
velocities  at  the  pressure  side  and  higher  velocities  at  the  suction  side.  This  velocity 
difference  translates  into  a  pressure  difference  across  the  passage,  which  results  in  a  flow 

turning  angle  of  1 1 .78°  at  the  trailing  edge. 

Figure  3.4d  shows  the  blade  wakes.  The  velocity  gradient  observed  across  the 
passage  width  is  not  seen  anymore  at  plane  7,  behind  thetrailing  edge. 

The  tip  leakage  vortex  has  an  additional  effect  on  the  flow.  The  tip  vortex 
produces  a  displacement  effect  that  constricts  the  flow  towards  the  trailing  edge  slightly 
decreasing  the  deceleration  of  the  flow.  For  the  present  flow,  with  a  tuning  angle  of 
11.78°  and  LU=25m/s,  continuity  equation  predicts  a  flow  speed  (on  average)  at  the 

trailing  edge  of  17.62m/s  (computed  as  \  9°)  ’  where  65-1  °  is  the  angle 

between  the  leading  edge  line  and  the  verticalsee  figure  2.2).  This  calculation  gives  a 
deceleration  ratio  of  29.5%  consistent  with  the  measured  28%. 
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Table  3.3  presents  the  maximum  mean  velocity  deficit  with  and  without  the  grid 
computed  as  (1-  Uc/Ucfree).  Where  Ucfree  is  the  local  center  streamline  aligned  velocity 
component  outside  the  tip  leakage  vortex. 


Table  3.3 


JC/Ca 

With  turbulence  grid 

With  out  turbulence  grid 

0.48 

0.45 

0.45 

0.77 

0.63 

0.61 

0.98 

0.76 

0.72 

1.26 

0.77  H 

0.74 

These  values  are  plotted  on  figure  3.5.  As  it  can  be  seen  the  velocity  deficit 
increases  to  the  trailing  edge  and  the  increment  is  a  little  greater  with  the  grid  on. 

On  figure  3.6a,  b,  c  and  d  mean  cross  flow  vectors  (computed  as  -J{K2+Wc2)) 

have  been  plotted  at  the  same  stations  before  mentioned.  Since  the  mean  velocity  field 
distribution  is  similar  with  and  without  the  grid,  only  one  set  of  plots,  that  corresponding 
to  the  case  with  the  grid,  is  shown.  Cross  flow  vector  plots  are  presented  in  the  same  way 
mean  stream-wise  velocity  plots  were,  with  the  exception  that  no  interpolation  has  been 
performed,  in  other  words  these  are  the  actual  measured  points. 

The  presence  of  the  tip  leakage  vortex  can  be  visualized.  As  shown  in  these 
figures  the  region  under  the  influence  of  the  tp  leakage  vortex  increases  as  the  flow  goes 
downstream,  reaching  a  height  of  approximately^  =0.4ca. 

Figure  3.7a  and  b  shows  the  variation  of  the  mass  flow  deficit  with  respect  to  the 
potential  core  per  unit  density  per  unit  width  (displacement  thicknessof  the  tip  leakage 
vortex  and  the  wakes)  with  and  without  the  grid.  The  displacement  thickness  was 
calculated  as: 

U,S=](U,-U)dy 

y\ 

where  Ue  is  the  local  axial  velocity  component  of  the  potential  core  at  the  edge  of  the 
viscous  region  (near  the  blade  wake  U e  was  computed  as  the  average  between  the  edge 
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velocities  on  both  sides  of  the  wake),  U  is  the  axial  velocity  at  either  the  vortex  or  the 
wakes,  yi  was  set  at  0.00254m  (minimum  height  at  which  measurements  were  taken) 
from  the  lower  end- wall,  yi  was  set  at  0.127m  from  the  lower  end  wall.  Figure  3.7a 
shows  the  presence  of  the  tip  leakage  vortex  at  planes  4  and  5.  Figure  3.7b  shows  the  tip 
leakage  vortex  at  planes  6  and  7,  the  boundary  layer  at  plane  6,  on  the  blade  4  side  and 
both  blade  wakes  behind  the  trailing  edge  at  plane  7. 

These  pictures  show  what  was  not  so  visible  on  the  previous  pictures.  The  size  of 
the  viscous  region  and,  therefore,  its  influence  are  increased  with  grid  on. 

The  last  two  sets  of  plots  involve  mean  stream-wise  vorticity  (£2^)  and 


turbulence  kinetic  energy.  Q^was  computed  as 


VQ 

FI 


where  V  is  the  velocity  vector 


and  Qis  the  vorticity  vector.  The  vorticity  components  in  axial  tangential  reference 
system  were  calculated  by  taking  the  curl  of  the  mean  velocity  vector.  To  calculate  the 
derivative  in  the  V  direction  it  was  assumed  that  the  derivatives  of  the  mean  velocfy 
components  in  the  mean  flow  direction  were  zero.  To  get  the  other  two  components  finite 
difference  approximations  were  performed.  Turbulence  kinetic  energy  was  computed  as 

(»?+y;2+^) 

2 

Figures  3.8a,  b  and  c  show  mean  stream-wise  vorticity  normalized  on  U„s  and  ca, 
at  x/ca=0.48,  x/ca=0.11,  x/ca= 0.98  with  the  grid.  Figures  3.8d,  e  and  f  show  the  same 
quantity  at  the  same  planes  but  with  out  the  grid.  It  is  clear  that  the  grid  reduces  the 
intensity  of  the  vorticity.  It  also  can  be  seen  a  wide  Egion  of  weak  negative  vorticity 
topping  a  small  region  of  positive  strong  vorticity  underneath.  This  negative  vorticity  is 
attributed  to  the  influence  of  a  secondary  vortex  spread  out  of  the  blade  near  the  leading 
edge.  See  Muthanna  (1998)  for  details.  Figures  3.9a,  b,  c  and  d  present  contours  of  the 
normalized  on  U2cS  with  the  grid  and  figures  3.9e,  f,  g  and  h  present  contours  of  the  with 
out  the  grid.  The  tip  leakage  vortex  is  shown  as  region  of  high  turbulence  whose  intensity 
is  increased  as  a  consequence  of  the  turbulence-generating  grid. 

These  set  of  plots  also  suggest  that  there  seems  to  be  a  region  above  y=0.55ca  not 
influenced  by  the  tip  leakage  vortex  along  the  passage.  Figures  3.10a,  b  and  c  show  the 
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distribution  of  C//LU  at  planes  5,  6  and  7,  figures  3.11a,  b  and  c  and  3.12a,  b  and  c  show 
V/Uoos  and  W/U~s  respectively.  These  contours  show  that  at  y=0.55ca,  indicated  by  a 
black  line  across  the  plots,  V  component  is  just  2-3%  of  the  other  two  components.  These 
region  extends  beyond y=l.lca.  The  upper  limit  was  not  determined  since  measurements 
were  made  up  toy=0.91ca. 

According  to  what  has  been  presented,  the  turbulencogenerating  grid  affects  the 
flow  field  in  the  following  ways: 

-)  It  barely  affects  the  structure  of  the  mean  flow  field  away  from  he  end- 
walls. 

-)  It  increases  the  size  of  the  zone  under  the  influence  of  the  viscous 
region. 

-)  It  increases  the  intensity  of  the  turbulent  flow  field. 

-)  There  exists  a  region,  not  influenced  by  the  tip  leakage  vortex  where  the 
flow  can  be  considered  two-dimensional.  According  to  the  data  the 
two-dimensional  region  starts  at  y=0.55ca  and  extends  beyond 
y=0.91ca. 

3.1.3  2-Dimensional  region 

According  to  the  previous  subsection  the  flow  field  can  be  considered  as  two 
dimensional  above  0 . 5 5c a  from  the  lower  end-wall.  The  comparison  between 
measurements  and  predictions  will  take  place  in  this  region  for  two  reasons.  The  first  one 
is  that  the  evolution  of  the  turbulence  distorted  by  the  blades  can  be  relatively  easily 
examined  and  modeled  here.  The  second  one  is  that  a  CFD  solution  of  the  mean  flow 
field  (RANS  calculations)  was  already  available.  Since  RDT  needs  the  mean  velocity 
field  to  model  the  distortion  and  recalliig  that  measurements  were  taken  at  only  8  planes, 
mean  velocity  measurements  do  not  provide  the  necessary  smooth  and  almost  continuous 
data  required  by  this  theory,  whereas  CFD  does. 

Figure  3.13  shows  contours  of  mean  velocity  magnitude  aty=0.73ca  (4").  The 
observer  is  looking  at  the  picture  from  above  and  the  flow  runs  from  left  to  right.  In  this 
picture,  the  flow  is  being  accelerated  over  the  suction  side  and  then  decelerated  as  the 
flow  goes  downstream.  The  velocity  contours  in  the  deceleration  region  ac  almost 


27 


Chapter  3.  Measurements 


perpendicular  to  the  flow,  this  being  due  to  the  large  stagger  angle  of  the  cascade,  which 
induces  the  flow  to  slow  down  sooner  on  the  pressure  side  than  on  the  suction  side.  It  also 
can  be  seen  the  rapid  change  in  color  of  the  contour  plot  which  gives  an  idea  of  where  the 
distortion  process  that  the  turbulent  flow  suffers  is  most  rapid. 

2 

Figure  3.14a,b  and  c  show  turbulent  velocity  components  normalized  on  U 
The  axes  of  the  plot  are  axial  tangential  but  the  turbulence  intensities  has  been  compted 
using  the  flow  aligned  coordinate  system.  The  distribution  ofw/  and  v/  is  uniform  along 
the  passage  presenting  a  little  decay  of  around  15%  toward  the  trailing  edge;  the 
distribution  of  the  third  component  wj  shows  a  significant  decay  toward  the  trailing 

edge.  This  also  can  be  seen  in  figure  3.15  which  shows  the  variation  of  turbulent  stresses 
components  aligned  with  the  flow  along  the  middle  passage  streamline  (see  Table  3.4  for 
coordinates  of  this  streamline),  w'2  is  also  suppressed  over  a  region  along  the  blade 

surface  and  0.2 ca  perpendicular  to  it.  This  is  consistent  with  the  idea  that  the  turbulent 
fluctuations  are  suppressed  by  the  non-penetration  condition  at  the  blade  surface. 

The  non-penetration  condition  can  be  modeled  through  the  introduction  of  a 
velocity  potential  that  modifies  the  velocity  field  but  that  does  not  alter  the  turbulent 
vorticity  distribution.  According  to  previous  studies  on  this  subject,  see  Wilker  et  al 
(1996)  and  Aronson  et  al  (1997),  the  suppression  of  the  normal  to  the  wall  component 
extends  over  a  region  of  about  one  integral  length  scale  in  thickness  from  the  surface  of 
the  boundary,  whereas  the  suppression  of  the  other  two  components  is  only  important 
within  a  region  of  about  one  order  of  magnitude  thinner.  This  explains,  perhaps,  why  the 
influence  of  the  non-penetration  condition  is  not  observed  on  bothw^  and  vs . 

Figure  3.16a,  b  and  c  show  the  stream- wise  integral  length  scales  distribution  over 
the  passage.  The  values  are  computed  in  flow  aligned  coordinates  system  and  equations 
3.4  was  used. 

Again,  the  distribution  is  fairly  uniform  over  the  passage  except  for  the  length 
scale  related  to  w's  which  has  been  significantly  suppressed  in  the  near  the  wall  region 
confirming  the  idea  of  the  additional  distortion  due  to  the  non-penetration  condition  or 
blade  blocking. 


28 


Chapter  3.  Measurements 


Figures  3.17a,  b  and  c  show  the  three  time  auto-spectra  components  in  the  stream- 
wise  direction  GUUJ(due  to  tO,  Gm(due  to  vs%  (due  to  measured  along  the 

middle  of  the  passage  streamline  at  six  different  locations  and  normalized  on if  Table 
3.4  shows  the  coordinates  of  these  locations.  Plane  1  is  not  considered  since  there  is 
almost  no  difference  in  the  spectrum  between  plane  0  and  1 . 


Table  3.4-Coordinates  of  the  center  streamline  in  axial  tangential  reference  system 


Center  streamline  z  /  ca 

coord. 

0 

OO 

o 

1 

1.69 

2 

0 

0 

3 

0.27 

-0.49 

4 

0.48 

-0.82 

5 

0.77 

-1.24 

6 

0.98 

-1.5 

7 

1.26 

-1.9 

Pictures  3.17a  and  b  show  a  little  variation  of  the  spectra  in  the  low  frequency 
region  while  the  variation  in  the  high  frequency  zone  is  more  important,  this  is  consistent 
with  the  uniformity  observed  on  the  integral  length  scales  on  figures  3.16a  and  b  and  with 
the  15%  decay  of  the  Reynolds  stresses  observed  on  figures  3.14a  and  b.  Figure  3.17c 
shows  a  significant  variation  on  both  the  low  and  high  frequency  regions  of Gwws  which  is 
consistent  with  the  variations  experienced  by  both  the  normal  to  the  wall  stress 
component  and  by  its  related  integral  length  scale  (see  figures  3. Be  and  3.16c). 
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4.1  Rapid  Distortion  Theory  (RDT) 

The  characteristics  of  a  turbulent  flow  can  be  altered  by  imposing  a  distortion  on 
the  stream.  For  instance,  the  contraction  section  of  a  wind  tunnel  reduces  the  turbulence 
levels  before  a  flow  enters  a  wind  tunnel  test  section.  Basically  the  imposed  distortion 
produces  large-scale  variations  on  the  mean  velocity  components  what  introduces  strains 
over  a  large  region  of  fluid. 

The  predictive  component  of  the  present  work  is  based  on  the  application  of 
Batchelor  and  Proudman's  (1954)  theory,  which  uses  the  linearized  version  of  Cauchy's 
equation  to  predict  distortion  effects.  Cauchy's  equation  is  the  Lagrangian  form  of 
vorticity  transport  equation: 


(Qi(d)  =  ^-(Qj(a) 


where  a)  (a)  is  the  vorticity  of  a  fluid  element  at  the  beginning  of  the  distortion  (time?’) 


when  its  position  is  determine  by  the  position  vector  a  ,to,  (a)  is  the  vorticity  of  the  same 
fluid  particle  at  time  t  during  the  distortion  when  its  position  is  determined  by  the  vector 


x . 


Finally 


is  rate  of  change  in  the  vorticity.  This  term,  known  as  the  distortion 


tensor,  will  depend  on  the  mean  velocity  field  in  the  neighborhood  of  the  particle.  It 


represents  the  modification  suffered  by  a  fluid  particle  as  consequence  of  the 
superimposed  distortion. 

For  the  linearized  form  of  Cauchy's  equation  to  be  applied  two  assumptions  are 
necessary  to  be  made: 

a)  The  duration  of  the  distortion  is  so  small  that  the  influence  of  viscous 
dissipation  can  be  neglected  during  the  process. 

b)  The  contribution  to  the  relative  velocity  of  two  neighboring  particles 
from  the  turbulence  should  be  negligibly  small  compared  to  that  of  the 
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distortion.  That  is  to  say  that  the  influence  of  turbulence  upon  itself 
should  be  negligible  small. 

In  other  words,  during  the  distortion  the  displacement  history  of  every  fluid  particle  is 
only  determined  by  the  history  of  the  distortion. 

Mathematically  this  is  expressed  as  follows: 

t 

a)  ^e(t)dt «  0.5  xwj2 

t' 

where  eft)  is  the  rate  of  dissipation  due  to  viscosity  andw/  is  the  Reynolds  normal  stress 
in  the  flow  direction. 

This  condition  could  also  be  expressed  as  TD  «  T (/) ,  where  TD  is  the  distortion 
time  (time  for  the  distortion  to  take  place)  and  x(l)  is  the  lagrangian  time  scaleof 

turbulence  computed  as  .  Here  A^s  is  the  longitudinal  integral  length  scale  at 
infinity  and  u'^s  is  the  turbulent  velocity  at  infinity  in  the  flow  direction. 

b)  S  >>u  /A  ,  S  is  the  strain  rate  of  the  distortion.  This  condition  implies  that 

/  oos  U°°S  7 

the  effect  of  the  superimposed  distortion  is  much  stronger  than  the  influence  of 
turbulence  upon  itself. 

According  to  Kevlahan  and  Hunt  (1996)  and  following  Hunt  and  Carruthers 
(1990),  these  assumptions  are  valid  when  the  linear  approximation  of  the  equation  of 
motion  is  used  to  compute  the  velocity  field  statistics  not  the  vorticity  statistics.  This 
could  be  because  they  are  not  based  on  the  vorticity  field  but  on  the  velocity  field. 

An  important  issue  in  applying  this  theory  is  to  get  the  distortion  tensor.  Batchelor 
and  Proudman  (1954)  assumed,  during  the  derivation  of  the  theory,  that  the  distortion 
tensor  was  homogeneous  and  that  it  varied  only  with  respect  to  time  (no  boundaries 
model).  So,  if  the  turbulent  field  was  homogeneous  before  the  distortion  it  remained  such 
after  the  distortion,  which  is  a  major  simplification  when  using  Fourier  transforms. 
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For  this  flow  configuration,  where  the  flow  passes  through  a  cascade,  the 
turbulence  is  homogeneous  and  isotropic  before  the  distortion  but  does  not  remain  the 
same  after  it,  since  the  distortion  is  not  uniform  (see  picture  4.1).  The  non-uniform 
distortion  generates  preferred  directions,  which  introduce  velocity  gradients,  keeping  the 
turbulent  field  from  being  both  homogeneous  and  isotropic.  However,  according  to 
Goldsteinb  (1979)  simplification,  if  the  integral  length  scales  of  the  turbulence  are 
smaller  than  the  characteristic  scale  of  the  distortion,  then  no  modifications  on  the 
method  are  necessary,  except  for  computing  the  distortion  tensor  at  every  point. 

For  the  present  work  a  different  approach  to  get  the  distortion  tensor  from  that  by 
Batchelor  and  Proudman  (1954)  was  used.  Darwin-Lighthill  drift  function,  see  Goldstein 
(1979),  seemed  to  be  a  simpler  way  to  get  the  distortion  tensor. 

Since  Batchelor  and  Proudman's  (1954)  model  was  a  non-boundary  model,  the 
effect  due  to  the  presence  of  boundaries,  called  non-penetration  condition,  was  not 
considered.  According  to  Hunt  and  Graham  (1976)  and  Hunt  and  Graham  (1997),  the 
non-penetration  condition  can  be  modeled  as  a  velocity  potential  that  modifies  the 
turbulent  velocity  field  but  not  the  turbulent  vorticity  field.  The  nonpenetration  condition 
will  affect  mainly  the  waHnormal  velocity  component  and  weakly  the  other  two 
components.  Perot  and  Moin  (1995),  Walker  et  al.  (1996)  and  Aronson  et  al.  (1997) 
showed  that  the  wall-normal  velocity  component  is  suppressed  over  a  region  about  an 
integral  length  scale  long  from  the  wall,  whereas  the  other  two  components  are 
suppressed  over  a  region  an  order  of  magnitude  thinner. 

Since  this  work,  as  the  initial  step  of  a  major  project,  is  concerned  about  the 
general  evolution  of  turbulence  convected  through  a  cascade  the  non-penetration 
condition  was  not  modeled  and  it  was  left  for  a  future  study. 

4.1.1  RDT  Conditions  evaluation 

As  it  has  been  stated  at  the  beginning  of  this  chapter,  for  RDT  to  be  valid  the 
history  of  the  turbulence  has  to  be  determined  only  by  the  history  of  distortion.  This 
condition  is  mathematically  expressed  by: 
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t 

a) \e(t)dt «  0.5  xw’2 

t' 

b)  S»uis/Aus 


For  the  present  flow  none  of  these  conditions  were  met.  The  first  condition,  and 
assuming  steady  state,  becomes  sAt,  where  At  is  the  period  of  time  for  the  distortion  to 
take  place.  It  was  estimated  as  c/LU;  c  is  the  airfoil  chord  equal  to  0.254m;  U„s  equal  to 
25m/s.  e  is  the  rate  of  dissipation.  Its  absolute  value  is  Q00018,  normalized  on  Ujlc  (see 
section  3.1.1  for  calculations).  After  evaluation,  the  left  hand  side  of  the  inequality  gives 
0.21.  The  right  hand  side,  for  an  inflow  turbulent  velocity  component  aligned  with  the 
flow  direction  of  0.62,  gives  0.31. 

In  case  of  the  second  condition,  the  strain  rate  of  the  distortion  S  has  been 


computed  by  a  finite  difference  approach 


along  the  middle  passage  streamline. 


The  left  and  right  hand  sides  of  the  inequality  31  and  31  respectively,  at  the  lade  row 
leading  edge  line.  At  the  trailing  edge  the  left  and  right  hand  sides  are  12  and  32 
respectively. 

Obviously  this  flow  lies  outside  of  the  range  of  validity  of  the  theory.  However, 
predictions  will  be  performed  for  one  of  the  objectives  is  to  checkthe  aptness  of  RDT  to 
predict  this  kind  of  turbulent  flow,  in  other  words  we  are  trying  to  push  the  theory  to  its 
limits. 


4.1.2  Theoretical  basis 

One  of  the  significant  components  in  RDT  calculations  is  the  computation  of  the 
distortion  tensor,  which  in  this  case,  represents  the  distortion  that  a  fluid  particle  suffers 
as  it  goes  through  the  cascade.  To  compute  the  distortion  tensor  Goldstein's  (1979)  high 
frequency  approach  was  chosen.  Since  the  integral  length  scales  of  turbulence  are  about 
20%  of  the  axial  chord,  Goldstein's  (1979)  approach  should  be  valid  except  in  the  vicinity 
of  the  blades,  where  blade  blocking,  not  accounted  for  in  this  approach,  dominates  he 
distortion  process. 
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Goldstein's  (1979)  model  used  the  drift  function  defined  by  DarwiHLighthill  as: 


_  *.  +  ] 

1 

1 

u„s  u„s  i 

uu(Xi,x2) 

w, 


(4.1) 


here  U\s  is  the  mean  velocity  parallel  to  the  flow  direction  at  location  1,-°°  represents  a 
non  distorted  point  sufficiently  away  upstream  of  the  distortion,*!  is  the  particle  actual 
position  vector  component  parallel  to  IU  and  X2  is  the  stream  function  normalized  on  the 
w 

free  stream  velocity: - . 

^  ooS 

Equation  (4.1)  represents  the  time  it  takes  a  particle  to  go  from  one  point  far 
upstream  on  a  given  streamline  to  another  point  on  the  same  streamline  In  other  words 
Xi  is  related  to  the  time  line  (see  Figure  4.2a). 

At  this  point  it  is  necessary  to  determine  the  relation  befcveen  Xi  and  X2  and 
Batchelor  &  Proudman's  (1954)  distortion  tensor.  They  defined  the  distortion  tensor  as: 


dt 


where  a7  is  the  initial  position  of  the  particle  at  timer',  *,  is  the  position  at  time  t,  Ui  is  the 
mean  velocity  componeit  and  Sy  is  the  Kronecker  delta.  This  equation  says  that  a 
particle  of  sides  (a„  aj)  at  initial  time  t’  will  change  its  shape  during  the  distortion 
according  to  Sy,  which  is  thus  the  distortion  tensor  represents  the  straining  rate  of  a  fluid 
particle  along  its  path. 

With  this  in  mind  if  we  think  of Xi  and  X2  as  determining  the  sides  of  a  particle  at 
t’  (see  figure  4.2b)  then 


’dxx 

dx{ 

ds 

dn 

dx2 

dx2 

.  ds 

dn  . 
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is  the  inverse  of  the  distortion  tensor.  Where  s  and  n  form  a  curvilinear  coordinate 
system,  s  along  the  streamline,  positive  downstream,  and/?  perpendicular  to  it.  So 


dxl 

dx{ 

ds 

dn 

dx2 

dx2 

.  ds 

dn 

(4.6) 


is  the  distortion  tensor. 

4.1.3  Evaluation  of  the  distortion  tensor 

4.1.3.1  Computation  of  the  derivatives 

S,2  component  is  the  hardest  to  compute  since  numerical  integration  along  a 
streamline  and  differentiation  in  n  has  to  be  performed.  Such  a  method  is  difficult  to 
accurately  implement  because  of  the  cumulative  error  in  integration. 

To  improve  the  accuracy  of  the  numerical  calculaions  the  problem  may  be  recast 
as  follows: 


hi  h±=\I-  .ds  (4.2) 

^  { us2 

Equation  4.2  represents  the  change  of  the  Drift  function  between  two  points  .4  and 
B  (see  figure  4.2a).  In  this  equation,  V  is  the  mean  velocity  vector,  Us  is  the  mean 
velocity  magnitude,  U^s  is  the  mean  velocity  in  the  free  stream  and  ds  is  a  differential 
streamline  element. 

If  the  change  between  two  other  points,  D  and  C,  is  considered  (see  figure  4.2a): 

,=  fJl.^  (4.3) 

U„s  du. 
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now, 


which  is  the  line  integal  around  loop  ABCD  (see  picture  4.2a).  If  we  remember  thatXy  is 

X  X 

related  to  the  time  line  and  that  U„s  is  homogeneous,  then  — 11  and  —j —  will  cancel  out 

each  other  since  they  have  the  same  value.  Now,  If  Stoke  fc  theorem  is  applied,  then  4.6 
becomes: 


1 

II 

t-s 

t^j‘ 

1 

<3 

X 

■dA 

=  -f 

IjL+v(vuj)xv 

J  TJ  J 

ABCD  w  s  Ar 

u. 

J 

Ar 

Us  J 

■dA  (4.7) 


where  Ar  is  the  area  within  loop  the  ABCD.  If  the  separation  between  two  streamlines  is 
sufficiently  small  then  dA  =  8n8sk  where  8n  is  the  distance  between  two  adjacent 
streamlines  and  k  is  the  area  unit  vector  perpendicular  to  n  and  s.  With  continuity 


equation  8n  =  8nA 


Therefore  equation  4.7  can  be  recast  into: 


_£iiL  =  _f[v(l luAxvXssSn^k  (4.8) 

TJ  U  J  1  '  U 

U~s  U°°s  A  s 
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assuming  irrotational  mean  velocity  field. 
Now, 


with  equation  4.8 


'a X ,  ^ 


dn 


)B 


1C  B 

Sn„ 


^  ^ n  )b 


*\C  B  _ 

8tln 


with  continuity  equation  8nA  =8nB  — —  equation  4.10  can  be  rewritten  as: 

oo.9 


=  u  UBs  f  v(l / K2  )x  —  •  kds 
dn  °°s  Bs{  V  F 


which  is  the  expression  to  compute  Sn  at  B. 
In  the  streamline  direction, 


d(X{/U„s)  JAU 


v  7 


ds 


ds 


since  V 'ds- Usds  : 


B  1 

af— ds 

8(X,/[/^)  JjJs  _  1 

8s  8^ 


&  */, 


(4.9) 


(4.10) 


(4.11) 


(4.12) 
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The  computation  of  the  derivatives  of  X2  is  simpler.  Recalling  that  X2  is  the 
streamfunction  normalized  on  the  free  stream  velocity  and  that  the  streamfunction 
represents  the  volumetric  flow  between  two  adjacent  streamlines,  then  AX2  =  U,SnlUm. 
This  formula  givesthe  value  of  X2  at  the  desired  point.  Then, 

_&Xj.  _  _  Us  (4.11) 

dn  5n  5n  U„s 


finally,  sinceXi  is  constant  along  the  streamline 


ds 


(4.12) 


which  gives,  finally  the  last  of  the  distortion  tensor  components.  Then 


dn 

Us 


(4.13) 


4.1. 3.2  R.A.N.S.  Solution  of  the  velocity  field 

The  velocity  field,  necessary  to  estimate  the  distortion  tensor,  was  calculated  by 
performing  two-dimensional RANS  calculations,  see  Shine/  al.  (1999),  on  a  high  density 
unstructured  grid.  This  kind  of  grid  has  the  advantage  of  its  flexibility  in  describing 
complex  geometries  such  as  a  highly  staggered  cascade.  A  SpalartAllmaras  eddy 
viscosity  model  was  implemented. 

Figure  4.3  shows  die  high  density  unstructured  grid  used  to  get  the  flow  field 
solution.  The  RANS  calculations  give  a  fairly  accurate  velocity  field  solution,  away  from 
the  lower  end- wall  region,  except  for  that  it  does  not  account  for  the  little  reduction  in  the 
deceleration  of  the  flow  due  to  the  flow  constriction  introduced  by  the  tip  leakage  vortex. 
The  flow  in  the  cascade  tunnel  experiences  a  speed  reduction  of  28%,  whereas  the  CFD 
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predicts  30%.  Given  this  level  of  accuracy,  CFD  was  chcsen  over  the  measurements 
because  it  provides  an  almost  continuous  solution  of  the  velocity  field,  which  was  not 
possible  otherwise  since  measurements  were  taken  at  only  eight  different  planes  along  the 
passage. 

4.2  Evaluation  of  the  spectra,  Reynolds  stresses  and  length  scales 

In  order  to  perform  all  these  calculations  a  Matlab  computer  code  was  written. 
This  code  is  divided  into  3  parts;  the  first  one  takes  the  flow  field  RANS  solution  and 
computes  the  distortion  tensor  at  the  desired  points.  In  order  to  apply  the  method  just 
described  a  total  of  15  streamlines  were  determined  for  the  present  flow.  If  the  desired 
point  was  not  on  any  of  these  streamlines  a  linear  interpolation  was  performed.  The 
second  part  calculates  both  the  undistorted  and  the  distorted  energy  spectrum  tensor. 
Finally  the  third  part  calculates  onedimensional  energy  spectra,  Reynolds  stresses  and 
Integral  length  scales. 

To  compute  the  undistorted  three-dimensional  energy  spectrum  tensor  the  von 
Karman  interpolation  formula  (Hinze,  1975)  was  used: 

*„(*„*,  - * — TB-k*.  -*’«„)  (414a> 

*2.* ,)  =  -C“T- - - — -j (4-14b> 

*“  -  (l +(klkjf 

- ! — (4-14c) 

k- k-  (i +(k/kjy 

where  O  is  the  three-dimensional  undistorted  energy  spectrum  tensor,  wMJis  the  turbulent 

velocity  component  in  the  flow  direction  at  infinity  (in  the  present  study,  infinity  is 
represented  by  the  flow  measured  at  the  first  plane),  ki,  k2  and  k$  are  the  wave  number 
vector  components,  k  =  ■yfk~kj  is  the  modulus  of  the  wave  number  vector, 
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ka  =0.75/ and  AUOOJ is  the  longitudinal  integral  length  scale  at  infinity^,  is  the 
Kronecker  delta.  Finally  C  is  a  constant  defined  as  follows: 


55  T(5/6) 

367r  2TVrcT(l/3) 


=  0.1156071356475;  f(5/6)  and  T(l/3)are  1.1288  and 


2.6789  respectively. 

These  formulae  will  be  used  to  describe  the  nondistorted  turbulence.  In  this 
study,  the  non-distorted  flow  field  is  described  by  the  data  measured  at  the  first 
measurement  plane.  Note  that  it  is  not  necessary  to  solve  for  every  component  of  the 
energy  tensor  but  only  for  the  main  diagonal  components. 


Equation  4.14  was  normalized  as 


follows: 


<D 


O, 


■ki, 


u„ 


<j)'  =  kl ,  =%^-k]s  giving  the  final  non-dimensional  expression  for  the  non 

u'L  M’W  U'L 

distorted  energy  spectrum  tensor: 


©,;«  x.o-  — -  k’2s, , ) 

(l+r1)* 

(4.15a) 

KXXX) - 

(l+r2)* 

(4.15b) 

(i +r!)‘ 

(4.15c) 

The  distorted  energy  spectrum  tensor  was  computed  by  applying  Batchelor  & 
Proudman’s  (1954)  formula: 


®l-disXX'W  =  £,U£mnpe  jabCuwlXl  (416) 


where  /= 1, 2,  3  and  j=  1, 2, 3. 


40 


Chapter  4.  Predictions 


X*,  =k](S*H)~x  is  the  non-dimensional  distorted  wave  number  vector,  Sy  is  the 

distortion  tensor,  is  the  unit  alternating  tensor  and  <!>*.(&*)  is  the  non-dimensional 

undistorted  energy  spectrum  tensor  computed  with  equations  4.15a,  b  and  c.  Note  that 
equation  4.16  can  solve  for  all  the  nine  components  of  the  tensor  but  only  the  main 
diagonal  components,  those  related  to  us',  vs'  and  ws',  will  be  used. 

Evaluating  one-dimensional  energy  spectra,  Reynolds  stresses  and  integral  length 
scales  from  the  distorted  three-dimensional  spectrum  function  involves  Inverse  Fast 
Fourier  Transforms  (IFFT).  The  full  three-dimensional  correlation  function  is  defined  as: 


Qy (x] , *2 , x3 , 0  =  J  J  J dx^x2dx^ihdist (XvXi’Xi, t)e[l 


i(Z\xi+X2x2+Xixi  )1 


(4.17) 


Equation  4.17  is  the  inverse  Fourier  transform  of  the  three-dimensional  energy 
spectrum  function.  Here  xi,  x2,  and  x2  define  a  certain  distance  from  the  point  where 
turbulence  is  being  considered  and  Oy_dlsl  is  the  three-dimensional  if  component  of  the 

distorted  energy  spectrum  tensor.  To  estimate  the  one-dimensional  spectra,  Reynolds 
stresses  and  integral  length  scales  only  the  zero  separation  parts  of  the  spectrum  are 
needed.  This  means  that  xj_  x2,  and  *5  are  zero,  which  implies  that  the  exponential  term  of 
equation  4.17  vanishes.  This  leads  us  to: 


«r (Xi>y\’z\)  =  j  (x')dxldxldx  1  (4. 1 8a) 

R  . 

x 

vr(Wi.*,)  =  J  ®22-<ns,(x’)dX}dxldx'l  (4.18b) 

R  , 
x 

wf(x,,y,,z,)  =  J  $*33-^ (x')dx\dx\dx\  (4.18c) 

R  * 


Equations  4.18a,  4.18b  and  4.18c  give  the  Reynolds  stresses  normalized  on  the 
inflow  turbulence  levels,  where  ,  d>22 _dist  and  033_dis,  are  the  main  diagonal 
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components  of  the  distorted  energy  spectrum  tensor  calculated  with  equation  4.16  and 
Rx.  is  the  non-dimensional  distorted  wave-number  space  formed  by  >  Xi  and  Xi 

Similarly,  the  one-dimensional  energy  spectrum  function  is  calculated  as: 

E  uus(X\)  =  2  j  j^uu-disl^X 3^X2  (4-19) 

-Xi  -Xi 


where  Euus*  is  the  one-dimensional  energy  spectrum  normalized  on  2  /  kes .  Similarly 
for  Eyyg  and  R u!w,v  • 

Finally,  once  the  one-dimensional  energy  spectrum  has  been  obtained,  the  integral 
length  scales  were  calculated  as: 


e^ksqtujx 


4  ul 

El 

,(2n50/Us)U„s 

4  ul 

E'^VtcSO/U'W,,, 

(4.20a) 

(4.20b) 

(4.20c) 


where  Kus,  A’vJand  Kws  are  the  integral  length  scales  normalized  on  kes  and, 
E’uus{2n5QIUs),  E'ms(2n50IUs)  and  £^(2&50/t/,)  are  the  spectral  energy  values  at 
the  stream-wise  wave  number  for  50Hz.  It  was  chosen  50Hz  instead  of  0Hz  for 
consistency  with  measurements. 

Performing  these  transforms  with  adequate  frequency  range  and  resolution  isat 
first  sight  a  near  impossible  computational  task  in  terms  of  memory  and  computational 
time.  This  situation  can  be  highly  improved  by  realizing,  first  that  only  the  zero 
separation  parts  of  the  spectra  are  needed.  This  implies  that  the  exponential  term 
vanishes.  Hence  the  inverse  fast  Fourier  transforms  are  reduced  to  simple  integrals  over 
the  wave  number  space  (see  equations  4.18  and  4.19).  Therefore  the  trapezium  rule, 
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instead  of  inverse  fast  Fourier  transforms,  can  be  applied  reducing  the  computationi 
time. 

Further  improvements  in  computational  time  and  memory  can  be  introduced  by 
using  semi-logarithmic  frequency  grid.  A  comparison  between  the  exact  solutions  of  the 
one-dimensional  energy  spectrum,  see  Hinze  (1975), 


E'{k\)  = 


'2*0.75 

7T 


(i 


+*,•7 


~5 


(4.21) 


and  the  one-dimensional  energy  spectrum  computed  by  means  of  the  trapezium  rule 
determined  that  the  necessary  wave  number  range  and  resolution  to  get  an  accurate 
answer  in  the  low  and  intermedate  frequency  range  was  ki  =±400  and  each  direction  of 
the  wave  number  space  divided  into  64  points.  A  grater  range  was  not  necessary  since 
von  Karman  formula  cannot  model  the  highest  frequency  region  of  the  spectrum.  Plus, 
measurements  in  this  region  wee  dominated  by  instrument  noise. 


4.3  Results  and  discussions  of  predictions 


As  mentioned  in  the  previous  section,  all  the  predicted  quantities  were  normalized 
on  turbulent  values  characterizing  the  incoming  flow  such  asw^ 2  and/or  A^ .  So  to  get 
the  values  shown  in  the  following  figures,  the  predicted  quantities  had  to  be  multiplied  by 
u^2  (in  the  case  of  Reynolds  stresses)  or  by  A^in  the  case  of  the  integral  length 
scales).  These  inflow  values  were  curvefitted  so  as  to  get  the  input  quantities  for  every 
streamline. 

Figure  4.4a,  b  and  c  show  the  distribution  of  RDT  predicted  Reynolds  stresses 
normalized  on  U2ccS.  As  it  can  be  seen  that  the  values  of  both  u  l2s  and  v*2*  increase  along 
the  passage  (especially  u’2 s),  which  is  totally  the  opposite  of  what  the  measurements 
illustrate  (see  figures  3.14a  and  b). 

wi2s  decreases  a  little  but  not  as  dramatically  as  the  measurements  show  where  the 
decay  is  of  about  50%  toward  the  leading  edge  and  near  the  blades  (see  figure  3.14c). 
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Figure  4.4c  does  not  show  either  the  suppression  w*s  suffers  near  the  blade,  however  this 
was  expected  since  blade  blocking  was  not  modeled.  Nevertheless  RDT  seems  to  account 
for  part  of  the  decay  experienced  by  the  v/,  along  the  passage. 

When  decay  due  to  viscous  effects  is  taken  into  account  in  the  calculations, 
predicted  values  get  better.  Viscous  effects  were  computed  in  a  very  sinple  way.  A  decay 
law  was  estimated  by  plotting  tke  values  measured  at  the  three  first  measurement  planes 
against*,; 

tkelU „  =0.001479e"°5819x> 

since  t=xs/U„s  at  the  inflow  region  (where  this  law  was  estimated),  then,  the  previous 
formula  can  be  rewritten  in  terms  of  timet. 

tke/U^  =  0.001479e-14-547' 


where  t  is  the  time  it  takes  the  particle  to  go  from  a  point  on  the  ^ero  time  line  plane  to  a 
different  point  downstream  on  the  same  streamline.  The  fcero  time  line’plane  was  placed 
at  *,=-0.1 6m,  between  the  first  and  the  second  measurement  planes,  perpendicular  to  the 
flow.  It  was  chosen  this  position  since  it  was  the  closest  to  £>lane  O’,  where  viscous  decay 
was  assumed  to  start.  The  first  measurement  plane  was  not  chosen  to  measure  timet  since 
a  plane  perpendicular  to  the  flow  was  needed,  otherwise  some  flow  particles  would  have 
been  more  viscous-effect-penalized  than  others. 

Figure  4.5a,  b  and  c  show  the  variation  of  Reynolds  stresses  due  only  to  viscous 
effects,  no  distortion  involved.  All  three  plots  present  exactly  the  same  pattern,  which 
was  expected  since  the  incoming  turbulent  flow  is  isotropic  and  viscous  decay  is  the  same 
for  the  three  components. 

Figure  4.6a,  b  and  c  show  Reynolds  stresses  contour  plots  where  both  viscous  and 
distortion  effects  were  combined.  The  distribution  of  predicted  values  ofw  ,  and  v  , 
(figures  4.6a  and  b)  is  quite  uniform  presenting  a  little  decay  toward  the  trailing  edge, 
especially  v'2s.  This  uniformity,  also  observed  in  measurements,  can  be  explained  since, 
according  to  these  calculations,  viscous  and  distortion  effects  seem  to  cancel  out  each 
other. 
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w^s  (figure  4.6c)  presents  a  stronger  decay  similar  to  that  observed  in 
measurements.  However,  as  the  non-penetration  condition  has  not  been  modeled,  w  s  is 
still  higher  than  it  should  be. 

Figure  4.7a,  b  and  c  show  the  predicted  values  of  integral  length  scales.  If  these 
quantities  are  compared  to  measurements  (see  figure  3.16)  it  can  be  seen  that  they  do  not 
agree  well.  This  difference  could  be  due  to  the  way  the  integral  length  scales  were 
computed.  They  are  estimated  as  a  function  of  the  lowest  frequency  bin  of  the  spectrum, 
where  Goldstein's  (1979)  high  frequency  approach  is  not  accurate.  This  approach  lies  on 
the  assumption  that  the  integral  length  scales  are  small  compared  to  the  main 
characteristic  length  (in  this  case  the  blade  chord).  The  integral  length  scales  are 
computed  in  terms  of  the  energy  spectra  values  for  the  lowest  frequency  bin  what  gives 
an  integral  length  scale  of  around  10%  of  the  blade  chord.  If  higher  frequency  binsare 
used  the  length  scales  will  be  smaller.  This  suggests  that  the  high  frequency  region  of  the 
spectra  is  in  better  shape  to  meet  the  condition  for  this  approach  to  be  valid. 

The  predicted  quantities  seem  to  be  smaller,  though  Aus  and  Aws  follow  the 
general  trend  observed  on  measurements. 

The  most  significant  differences  are  observed  on  Avs  where  RDT  predicts  a 
notorious  decrement  as  the  flow  reaches  the  trailing  edge.  Another  difference,  expected 
though,  is  seen  in  Aws  in  the  near-the-wall  region  due  to  the  non-consideration  of  non¬ 
penetration  condition. 

Reynolds  stresses  as  well  as  Integral  length  scales  plots  present  linear  features, 
this  been  due  to  the  way  the  input  data  was  handle  to  compute  the  turbulent  quantities 
distribution  along  the  passage. 

Figure  4.8a,  b  and  c  show  the  measured  and  the  predicted  spectra  (displaced 
100Hz  to  the  right).  These  values  have  been  normalized  on  LU  and  Au~s.  The  predicted 

quantities  of  the  spectra  were  normalized  on  ’  In  order  t0  get  the  new 

u'2 

normalization,  predicted  spectra  were  multiplied  by^2  ^  —  •  Note  that  the  local  value  of 

the  flow  aligned  Reynolds  stress  was  used.  This  way  the  viscous  effects  were  roughly 
included  in  the  spectra.  Predicted  quantities  have  been  shifted  to  the  right  by  multiplying 
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UoosAuoos  times  10.  As  it  can  be  seen  predictions  reproduced  the  general  trend  observed  in 
measurements. 

According  to  what  has  been  presented  Rapid  Distortion  Theory  has  been  able  to 
reproduce  the  gross  features  of  this  flow.  The  results  achieved  are  encouraging  if  we 
consider  the  fact  that  the  present  flow  lies  totally  outside  of  the  range  of  validity  of  the 
theory. 
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Measurements  of  a  turbulent  flow  field  past  a  linear  compressor  cascade  were 
made  by  using  hot  wire  anemometry.  Basic  RDT  predictions  were  performed  in  an 
attempt  to  determine  the  range  of  application  of  this  theory  when  used  to  predict  the 
effects  of  distortion  in  a  two  dimensional  region  in  a  highly  staggered  cascade  flow. 

Based  on  the  results  from  measurements  and  RDT  predictions  the  following 
conclusions  can  be  made: 

•  Measurement  reveals  the  formation  and  evolution  of  the  tip 
leakage  vortex.  The  tip  leakage  vortex  is  represented  by  a  larger 
region  of  positive  stream-wise  vorticity  bounded  on  top  by  a 
smaller  region  of  weak  negative  stream-wise  vorticity. 

•  The  turbulence  grid  increases  the  turbulence  levels  and  it  increases 
the  size  under  the  influence  of  the  viscous  region. 

•  In  the  two-dimensional  region  us'2  and  present  a  uniform 
distribution.  RDT  calculations  suggest  that  the  distortion  suffered 
by  the  flow  increases  these  quantities  in  almost  the  same  amount 
that  they  decrease  due  to  viscous  dissipation. 

•  In  the  case  of  ws'2,  viscous  and  distortion  effects  combine  together 
so  that  ws'2  decreases  by  about  50%  towards  the  trailing  edge  line. 
ws'2  and  its  related  integral  length  scale  experience  a  strong 
suppression  towards  the  blade  walls  due  to  the  non-penetration 
condition  effect. 

•  RDT  calculations  plus  viscous  effects  considerations  were  able  to 
provide  the  general  trend  of  the  evolution  of  the  turbulent 
fluctuations  through  the  cascade,  especially  of  us  ,  vs  and  ws 
components.  In  this  case  predictions  were  within  10%  of  the 
measured  values. 

•  However,  important  differences  can  be  observed,  vv,'2,  in  the  near- 
the-wall  region,  experienced  an  important  suppression  due  to  non¬ 
penetration  condition,  not  observed  in  predictions.  The  same  thing 
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happens  to  the  related  lateral  length  scale  AW|S.  These  differences 
were  expected  since  blade-blocking  effects  were  not  modeled. 

•  The  most  important  differences  can  be  seen  in  the  integral  length 
scales.  For  Aus  and  Aws  predictions  provide  the  general  trend  but 
the  magnitudes  are  within  20%.  For  Avs,  RDT  predicts  a  dramatic 
decay  through  the  passage  (45%-50%),  not  observed  in 
measurements  (the  decay  is  about  1 5%-20%).  These  differences 
may  rest  on  the  way  integral  length  scales  were  computed. 

The  agreement,  though,  is  more  than  expected  if  the  fact  that  the  analyzed  flow 
lies  totally  outside  of  the  range  of  validity  of  the  theory  is  considered.  In  other  words,  the 
this  study  suggests  that  a  proper  RDT  model,  including  viscosity  models  and  blade 
blocking  effects  could  predict  with  a  good  level  of  accuracy  the  main  characteristics  of  a 
turbulent  flow. 
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Fig.  1.1:  Turbulent  flow  being 
distorted  by  a  Rotor-stator  system 


Sidewall 


Exit  screens 


Fig  2.2:  Top  view  of  the  Virginia  Tech  Cascade 
Wind  Tunnel  Test  Section.  Dimensions  in  meters 


Traverse 

System 


Fig  2.4:  Cascade  Tunnel 
inlet  section  roof  top 


Fig  2.6:  Screens  at  the  exit 
plane  of  the  tunnel. 


Fig  2.7;Schematics  of  the 
aluminum  pressure  tabs  (not 
on  scale). 
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Fig  2.9:  General  View  of  the 
Cascade.  Passage  4,  where 
measurements  have  been 
taken  is  indicated  by  the 
white  arrow 
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Adjustable 


Fig  2.10:  Cross  section 
of  the  suction  slots 
(not  on  scale) 
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Fig  2.11:  Bleeding 


Ceramic  Tubes 
(length  -  10mm) 


Probe  Prongs  Stainless  Steel  Casing 

(length  -  40mm)  (length  =  63mm,  diam.  -  4.3mm) 


Delrin  Block,  (length  =  50mm, 
7.9mm  square  cross-section) 

ft  \ 


DUPONT  BergCon 
Type  Connector 


Electrical  Leads 
(length  *-  50mm) 


Fig  2.12:  Four  Sensors  Hot-Wire 
Probe 
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Fig  2.13:  Hot-Wire  Probe  in  the 
calibration  equipment 


Fig  2.15:  Probe 
holders 
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Fig  2.18:  Final  pressure  distribution  after  calibrating 
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Fig.  3.1:  Scheme  showing  a  blade  passage,  the  8 
measurement  planes,  the  axial-tangential  and  center 
stream  line  aligned  reference  systems 


Endwall  boundary  layer 


Fig.  3.2:  Scheme  of  he  viscous  region.  The  blade  wakes, 
lower  end-wall  boundary  layer  and  tip  leakage  vortex 
can  be  visualized. 
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Fig  3.3a:  Ratio  between  turbulent  velocity 
components  at  0.0762m  off  the  lower  end- 
wall 
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Fig  3.3b:  Ratio  between  turbulent  velocity 
components  at  0.1016m  off  the  lower  end-wall 
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tt/v  /-  u/w  /- distribution 


Fig  3.3c:  Ratio  between  turbulent  velocity 
components  at  0.127m  off  the  lower  end-wall 


tkeAJ,J -distribution 


z/ca 

Fig  3.3d:  Averaged  turbulent  kinetic 
energy,  normalized  on  plotted 
along  the  first  measurement  plane 
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Fig  3.3e:  Comparison  between  von  Karman  interpolation 
formula  and  measurements  at  the  inflow  plane  of  the  three 
autospectra  components.  The  case  with  the  grid  is  presented. 

The  spectra  were  normalized  on  UKS2 
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Fig  3.3f:  Averaged  longitudinal  and  lateral  integral  length  scales  (in 
meters)  plotted  along  the  first  measurement  plane  with  the  grid  on. 
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e)  x/ca= 0.48 


IllllllllllllilllillU 


Fig.3.6b:  Cross  flow  vector  plot  at  x/c  =0.77  normalized  on 


Fig.3.6c:  Cross  flow  vector  plot  at  x/ca=0.9S 
normalized  on 


U 


z./c„ 


Fig.3.6d:  Cross  flow  vector  plot  atx/cfl=1.26 
normalized  on 


Fig.3.7a  Fig.3.7b 

Fig.3.7:Displacement  thickness  measured  a  the  last  four  measurement  planes 
with  and  without  the  turbulence-generating  grid.  Fig  4.7a  shows  the 
displacement  thickness  at  jc/ca=0.48  and  at  x/ca-d.ll  and  Fig.  4.7b  shows  the 
displacement  thickness  at  x/ca=0.98  and  at  x/c=  1.26 


74 


1 

0.8 

^0.6 
^  0.4 
0.2 
0 


e)  x/ca=OA8 


J&im 


inliiiiliiiiliiiiliiiiliniliiiillillllllllllllllllTI 
-2.8 -2.6 -2.4 -2.2  -2  -1.8 -1.6 -1.4 -1.2  -1  -0.8 -0.6 -0.4 -0.2  0 


z/c 


f)  x/c= 0.77 


l 

0.8 
o0.6 
^  0.4 
0.2 
0 


inn  mu  nm  mn  iiivfiTiiin  inn 


1 1 1 1 1 1 1 1 1 1 1 1 1 


-2.8 -2.6 -2.4 -2.2  -2  -1.8 -1.6 -1.4 -1.2  -1  -0.8  -0.6  -0.4  -0.2  0 
z/c 


0.019 
0.019 
0.018 
0.01? 
0.017 
0.016 
0,01 5 
0.015 
0.014 
0.013 
0.013 

10.012 
0.011 
0.011 
0.010 
0.009 
0.009 
0.008 
0.007 
0.007 
0.006 
0.005 
0.005 
0.004 
0.003 
0.003 
0.002 
0.001 
0.001 
0.000 


g)  x/c= 0.98 
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h)  x/c=  1.26 


Fig.  3.9:  Contour  plots  of  the  normalized  on  Uxs2 
with  out  the  grid. 
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Fig.  3.10:  Contours  of  the  mean  velocity 
component  in  axial  direction  normalized  on 
at  three  different  x  locations 
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Fig.  3.12:  Contours  of  the  mean  velocity 
component  in  tangential  direction  normalized 
on  £LC.  at  three  different  x  locations 


Fig.  3.14b:  Contour  plot  of  v/ 
normalized  on  U^f, 
at  y/c.  =  0.73  with  the  turbulence  grid. 


Fig.  3.14c:  Contour  plot  of  w2  8 
normalized  on  U„s2, 
at  y/ ca  =  0.73  with  the  turbulence  7 
grid. 
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Fig.  3.15:  Variation  of  the  Reynolds 
stresses  along  the  middle  passage 
streamline. 

Values  normalized  on  U  2 


Fig.3.17b:  Gvvs  autospectra  component  normalized  on  Ums2,  j 
different  locations  along  the  middle  passage  stream-line 


Fig.3.17c:  Gwm  autospectra  component  normalized  on  £/„/, 
different  locations  along  the  middle  passage  stream-line 


Fig  4.1:  Fluid  particle  suffering  non- 
homogeneous  distortion  due  to  a  blade  row 
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Fig.  4.2  a:  Schematic  for  the  Drift  function 
calculation 


Fig.  4.2b:  Relationship  between  a  particle  sides, 
the  Drift  function  and  the  Streamline 
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Fig.  4.3:  Close  up  view  of  the  grid  at  the  leading  edge  of  blade 
4  for  the  RANS  calculation.  The  negative  sign  in  the  z/c  axis 
label  is  to  adapt  the  R.A.N.S  calculation  reference  system  to 
the  tunnel  reference  system. 
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Fig.  4.4a:  RDT  predictions. 

component  of  the  Reynolds  stresses, 
normalized  on  V 2OOS. 


Fig.  4.4b:  RDT  predictions. 

v>2s  component  of  the  Reynolds  stresses, 
normalized  on  U2 _ 


U.8 


Fig.  4.4c:  RDT  predictions. 

component  of  the  Reynolds  stresses, 
normalized  on  V 2„s. 


Fig.  4.5a:  Decay  of  u’2s  due  only  to 
viscous  effects  normalized  on  V2^.  (no 
distortion  considered) 


Fig.  4.5b:  Decay  of  v’2s  due  only  to 
viscous  effects  normalized  on  U 2oes.  (no 
distortion  considered) 
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Fig.  4.5c:  Decay  of  w’2s  due  only  to 
viscous  effects  normalized  on  V 2eos.  (no 
distortion  considered) 


Fig.  4.6a:  RDT  predictions. 

,2S  component  of  the  Reynolds  stresses, 
normalized  on  V 2„s. 

Viscous  decay  has  been  considered. 
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Fig.  4.6b:  RDT  predictions. 

v,2s  component  of  the  Reynolds  stresses 
normalized  on  U2^,  viscous  decay 
considered 


Fig.  4.6c:  RDT  predictions. 

s  component  of  the  Reynolds  stresses, 
normalized  on  V 2W,  viscous  decay 
considered 
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Fig.4.8a:  Comparison  between  measurements  and  predictions  (decay  included)- 
Euus  spectrum  component  normalized  on  A^U^  plotted  against  kj*  A UMS  The 
predicted  quantities  have  been  shifted  by  multiplying  k1  Auoos  by  10. 


Fig.4.8b:  Comparison  between  measurements  and  predictions  (decay  accounted 
for)-2?vra  spectrum  component  normalized  on  Auoo*UMS  plotted  against  k,*  Auoos. 
The  predicted  quantities  have  been  shifted  by  multiplying  kj  Auoos  by  10. 
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Fig.4.8c:  Comparison  between  measurements  and  predictions  (decay  accounted 
for)-^^  spectrum  component  normalized  on  Auoos  U„s  plotted  against  kj  Aucos. 
The  predicted  quantities  have  been  shifted  by  multiplying  kj  Aucos  by  10. 
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